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Preface 


This monograph deals with the problems of mathematical physics which 
are Improperly posed in the sense of Hadamard. 

The first part covers various approaches to the formulation of improperly 
posed problems. These approaches are illustrated by the example of the classical 
improperly posed Cauchy problem for the Laplace equation. 

The second part deals with a number of problems of analytic continuations 
of analytic and harmonic functions. 

The third part is concerned with the investigation of the so-called inverse 
problems for differential equations in which it is required to determine a dif- 
ferential equation from a certain family of its solutions. 


Novosibirsk June, 1967 M. M. LAVRENTIEV 
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Chapter | 


Formulation of some Improperly Posed Problems 
of Mathematical Physics 


§ 1. Improperly Posed Problems in Metric Spaces 


The notion of correctness* introduced at the beginning of our century by 
the French mathematician HADAMARD plays an important role in the investiga- 
tion of the problems of mathematical physics. One often says that a problem is 
solved if its correctness is established. 

Various authors present notions of correctness which coincide in their 
essence but differ in details. We give one of the possible definitions of correct- 
ness which is convenient for our aims. 

Let ®, F be some complete metric spaces, and let Ag be a function with the 
domain of definition ® and the range of values F. Consider the equation 


Let us point out that most problems of mathematical physics can be reduced 
to the investigation of the solution of equation (1.1) with a given function A 
and right-hand side f. We say that the problem of solving (1.1) is properly 
posed if the following conditions are satisfied: 

1) The solution of (1.1) exists for any fe F. 

2) The solution of (1.1) is unique in ®. 

3) The solution of (1.1) depends continuously on the right-hand side f. 


In other words, the problem of solving (1.1) is properly posed if there exists 
a function Bf defined and continuous over all of F, which is inverse to the 
function A@. 

Linear problems are most often considered in mathematical physics. In 
this case ®, F are BANACH spaces, and A is a linear operator. The BANACH 
spaces ®, F in concrete problems are the known functional spaces C!, L,, We 
HAH s Sp, --. with the carriers in some n-dimensional space of the independent 
variables or on any part of the spaces of independent variables. 


* Translator’s note: Hereafter we shall refer to the correctness or incorrectness of 
problems which are respectively properly posed or improperly posed. 
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The first requirement of correctness is that the problem should not be 
overdetermined, and superfluous conditions should not be imposed. 

The second requirement is that the solution be unique. 

The third requirement of correctness, continuity of the inverse function Bf, 
arises from the fact that in the real problems of mathematical physics the right- 
hand side of equation (1.1) is obtained from measurements made with the aid 
of actual instruments and is therefore known only approximately. Therefore, 
it has been felt for a long time that if at any point f the function Bf is dis- 
continuous, then the solution g cannot be uniquely recovered from the right- 
hand side f- 

HADAMARD introduced the notion of correctness by giving an example of 
an improperly posed problem which became classical and was included in 
most text-books on mathematical physics. The example is the CaucHy problem 
for the LAPLACE equation. It is well-known that the solution of this linear 
problem does not depend continuously on the data obtained from any of the 
functional spaces mentioned above. On the basis of this HADAMARD concluded 
that CaucHy’s problem for the LAPLACE equation and, in general, all problems 
exhibiting a similar dependence of the solution on the right-hand side, do not 
correspond to any real formulations, i.e., they are not problems of mathematical 
physics. 

It was discovered later that HADAMARD’s conclusion was erroneous and 
many real problems of mathematical physics lead to problems which are 
improperly posed in the sense of HADAMARD. In particular a number of im- 
portant problems of geophysics lead to the CAucHy problem for the LAPLACE 
equation. 

There are still further examples of linear as well as nonlinear improperly 
posed problems which are important in the applications, namely the solution 
of the heat equation for negative time and Caucuy data on the boundary, the 
nonhyperbolic Caucuy problem for the wave equation, inverse problems of 
potential, and a number of inverse problems for differential equations. 

At present there exists a number of approaches to the investigation of 
improperly posed problems. We will now explain them using the above classical 
Caucny problem for the LAPLACE equation as an illustration. 

We consider one of the simplest versions of the CAucuy problem for the 
LAPLACE equation. Let u (x, y) be a twice continuously differentiable function 
of the variables x, y in the rectangle 


O<x<z, O<y<H (D) 
satisfying the following conditions 
Au=0, (x,y)eD (1.2) 


0 
5p BOD aua =o: 
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It is required to determine its values on the segment y = h, 0 < x < a by 
its values on the segment y = 0,0 < x < x. 


The problem stated is equivalent to the solution of equation (1.1) where 
~, f stand for the functions of the variable x 


gp =u(x,h) 
f =u(x,0). 


A is a linear integral operator with the kernel 


Ky (x, g)= 


aly 


> cosh”! kh-sinkx:-sinké. 
h=1 


The functional spaces ©, F considered in the given example are the HILBERT 
functional spaces with the scalar products of the form 


(91.92)= 1 Agi G2 ped 


o (1.3) 
f= 2X mit Sa S;EF 


where oy and f* are the FOURIER coefficient of the functions g,; and f;, and A, 
and 4, > O are some sequences. 


We note that spaces of such form are, in particular, the spaces w! with the 
norm 








oe 
Ox y 








! 
wlla: = x C; 
j=0 


Lz 


It may be easily seen that the solution of (1.1) in the case under considera- 
tion is unique, but for its existence and continuity it is insufficient to let ® and 
F be some W3 with finite /. 

To make the problem properly posed one can introduce a sufficiently 
‘‘strong’’ norm in the data space F or a sufficiently ‘‘weak”’ one in the space ®. 
It is quite evident that the solution of (1.1) is correct in the sense of our defini- 
tion, if the sequences A,, 4, in (1.3) satisfy the inequality 


Ae <C ty, (1.4) 


where C is a Constant. 


The inequalities (1.4), in particular, are satisfied by the following sequences 
Ay, 4, Corresponding to the above possibilities 


A=e*, b=1 


5 
A=1, =e”. aa 
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However, the above approaches have the following defects. In the first 
case the statement does not cover the range of problems in which the errors 
in the data may be considered to be small only in the “ordinary” functional 
spaces; in the second case solutions of the problem may be objects which are 
not functions in the usual sense and do not correspond to physical reality. 

We now formulate some approaches to the question of correctness of 
problems of the type under consideration which are free of the above defects 
and, in our opinion, quite natural from the standpoint of applications. 

The first approach consists of changing the notion of correctness, namely, 
to one having requirements different from 1), 2), and 3). 

We now state the requirements. In addition to the spaces ® and F and the 
operator A, let there be given some closed set Mc ®. 

We call the problem for the solution of (1.1) properly posed according to 
Tykuonov if the following conditions are fulfilled. 


1) It is a-priori known that the solution @ exists for some class of data and 
belongs to the given set M, pe M. 
2) The solution is unique in a class of functions belonging to M. 


3) Arbitrarily small changes of the right-hand side of f which do not carry 
the solution @ out of M correspond to arbitrarily small changes in the solu- 
tiong. 

We denote by M, the image of M after the application to the space ® of the 
operator A. 

Requirement 3) can be restated in the following manner, 

3) The solution of equation (1.1) depends continuously on the right-hand 
side fon the set M,. 

If M is a compact set the following statement holds (see [17]). 

If equation (1.1) satisfies the requirements 1), 2) of correctness due to 
TYKHONOV, then there exists a function « (t) such that 

a) « (t) is a continuous nondecreasing function with « (0) = 0. 


b) for any @;,@2 € M satisfying the inequality 


@(Ag,,Ag2)<e 
the following holds 


2(91,2)<a(e). 


Thus, the requirement of continuous dependence 3) is satisfied if 1) and 2) are 
satisfied. 


We note that, if a problem is properly posed according to TycHONOV and 
we replace the metric spaces ®, F by their supspaces M, M 4 then the problem 
becomes properly posed in the usual sense. 
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The necessity of examining spaces ®, F together with M, M, is due to the 
fact that in real problems the errors committed in the determination of the 
right-hand side f usually lead to f outside of My. The consideration of the 
problem according to TYKHONOV’s formulation gives the possibility of construct- 
ing an approximate solution with a certain guaranteed degree of accuracy in 
spite of the fact that an exact solution of (1.1) with approximate data either 
does not exist at all or may strongly deviatc from the ‘‘true”’ solution. 

In the Caucny problem for the LAPLACE equation we consider as the sct Mf 
the sets defined in the following manner: we denote by A n, 2n integral operator 
with the kernel 


alnN 


of) 8 


K,, 0%, O= cosh”! kh, sinkx-sinké 


and let 


epeM(I,h,,C) 
if 
p=Ay,,, lly a<C, (hy =H—h). 


We now give evidence that our formulation is natural for the case in which 
the CAucny problem for the LAPLACE equation is used for the solution of the 
geophysical problem of interpreting the gravitational or magnetic anomalies. 
The problem of interpretation of geophysical data is as follows: the characteris- 
tics of some physical field interacting with the inner layers of the earth are 
measured on the earth’s surface. Certain characteristics of the structure of 
these layers are required to be determined. In interpreting the constant magnctic 
and gravitational fields there often arises the following situation. An upper 
layer of the earth’s crust consists of uniform sedimentary rocks. The thickness 
of the sedimentary rock layer may be estimated from below by some constant 
M. The anomalies observed on the earth surface, i.e., the CAUCHY data for a 
harmonic function, are generated by some bodies lying at a depth exceeding H 
and the ultimate aim of an observer is to determine these bodies. This problem 
is rather complicated; its solution is not unique if one proceeds only from 
measurements of a field on the surface. However, by proper evaluation of the 
general geological situation and with well chosen supplementary hypotheses, 
the observer often succeeds in finding solutions with satisfactory accuracy. 

It is important in this problem to make good estimates on the number of 
the bodies, their approximate position, and their size. For the case in which 
the average depth of embedding of the bodies is substantially less than their 
horizontal dimensions and the spacing between them, the distinct extrema of 
the anomaly correspond to distinct bodies. However, if the average dcpth of 
imbedding roughly coincides with or excceds the characteristics horizontal 
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dimensions, then there is no such correspondence. It is possible for two bodies 
to cause one extremum situated in the gap between them. In order to obtain 
the case in which the depth of embedding is less than the horizontal dimensions 
one uses the solution of CAUCHY’s problem for the LAPLACE equation. The 
existence of the solution is ensured by the fact that the field to be measured is 
generated by real objects and by the fact that in the region of evaluation, the 
depth of which does not exceed H, singularities of the field do not occur. The 
fact that the solution belongs to the set M is guaranteed by the fact that in the 
earth’s crust one does not encounter bodies of excessive density or with magnetic 
intensity exceeding some completely determined constants. For example, when 
u(x, y) is a gravitational potential the above considerations guarantee tliat 
the solution belongs to a set M (2, A,, C), in the case when u is a component 
of the gravitational field stress— (1, A,, C), and in the case when w is a 
component of the magnetic field stress— M (0, A,, C). 

This approach to the question of the formulation of improperly posed 
problems was originally proposed by TyxHoNov [17]. A systematic investiga- 
tion of various problems, properly posed in the sense of TyKHONOV, was carried 
out in [28], [31], [50], [56]. A detailed description of the definitions as well as 
the investigation of a number of problems which are properly posed in the 
sense of TYCHONOV are contained in the monograph [35]. 

Another possible approach to the study of improperly posed problems is 
to let the notion of correctness remain unchanged, but change the notion of a 
“solution” of the problem. Besides the spaces ® and F and the operator A let 
there be given, as in the preceding case, some set M. By a quasisolution of 
equation (1.1) we mean an element @ defined by the relation 


e(Ag, f)=min @ (4g, f). 


We note that the problem of finding a quasisolution is one of nonlinear 
programming. From the theorems of nonlinear programming one obtains the 
corresponding theorems on the correctness of the problem of finding a quasi- 
solution. 

The problems which may be regarded as properly posed according to 
TYKHONOV are those for which the notion of quasisolution may be suitably 
introduced. In some cases the quasisolution may be non-unique, but when the 
problem is properly posed due to TYKHONOV, any quasisolution constructed 
using an approximate right-hand side will be one of the best versions of an 
approximate solution of the problem. 

The notion of a quasisolution was introduced in [45]. Even though, as 
previously mentioned, the range of applicability of this notion coincides with 
that of the notion of correctness due to TYKHONOV, the notion differs by offering 
greater clarity and simplicity. 
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The most general approach to improperly posed problems was given in the 
recently published works of TyKHONov. Following [18] we define a regularizor 
for Eq. (1.1) to be any one-parameter family of operators B, having domain of 
definition F and range of values @ and satisfying the following conditions 

1) for any t > 0 the operator B, is defined and continuous over all of F. 

2) for anyge®@: 

lim B,Ag=o. 
t70 

We call the problem of solving (1.1) regularizable if there exists a regularizor 
for (1.1). It may be easily scen that the regularizor permits us to construct a 
solution with a guaranteed degree of accuracy according to the approximate 
right-hand side. In fact, let f be the right-hand side of (1.1) with an error e, 
1.€., 


of, f)<e. 
We denote 
6,=B.f; 9,=B,Ap 


and estimate the quantity 0 (¢,@,). 
By the triangle inequality 


0(¢, 6)<0(9,9,)+0(Gn G,)- (1.6) 
It is clear that 
2(P, 9 )S%,(T 
(@ 0) 1 (t) (1.7) 
2 (Ps GB.) SM (T, €) 


where «, (t) is a function characterizing the rate of the convergence of B,Ag 
to @, and «, (tT, €) is the modulus of continuity of the operator B, at the point 
Ag. The functions «, (rt) and «, (t, €) for fixed t are continuous and monotonic, 
and x, (0) = «, (1, 0) = 0. 

Substituting (1.7) into (1.6) and setting t equal to a root of the equation 


Gt, (tT) = a2 (T, £) (1.8) 
where &, (t), is the inverse of «, (t), we obtain 
o(¢, @,) S20, (t). (1.9) 


It is evident that the root of equation (1.8), and therefore, the right-hand 
side in the inequality (1.9) tends to zero as € > 0. 

In contrast to the first two approaches, in the notion of regularizability we 
consider no additional sets other than the basic spaces ®, F. This accounts for 
the greater generality of this approach. Apparently, all the problems of mathe- 
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matical physics connected with real phenomena are regularizable. (HADAMARD 
regarded all such problems to be properly posed.) In [22] the regularizability 
of a rather wide class of equations (1.1) has been proved. 

In conclusion we point out the relation between the notion of a quasisolu- 
tion and a regularizor. Suppose the solution of (1.1) is unique. Let M, be a 
one-parameter family of compact sets (rt > 0) satisfying the requirements 


M,,¢M,,, ™1<t,U M,=®, 
t>O0 
and let B, be a continuous operator with domain of definition F such that B, f 
is a quasisolution of (1.1) with the set 4/ = M,,. It is evident that the operator 
B, is a regularizor for cquation (1.1). 
Besides the above three approaches to improperly posed problems, an 
approach related to the theory of probability has been developed to some 
extent. The next section is devoted to this approach. 


§ 2. A Probability Approach to Improperly Posed Problems 


The theory of probability was first utilized in the investigation of improperly 
posed problems in a paper by KHALFIN and SUDAKOV [66]. The authors restricted 
their investigation to the classical CaucHy problem for the LAPLACE equation. 
The main content of the paper is as follows. 

The Caucnuy data are assumed to consist of a sum of two items, the “‘exact 
data” and the “‘error”’. The “error” is a realization of a stationary random 
process. The question to be investigated is, under what conditions on the 
autocorrelation function of the error will this error insignificantly affect (in the 
average) the solution of a problem. It turns out to be sufficient that the spectral 
function of the error should be small in the HILBERT space with the exponential 
weight. 

KHALFIN and SUDAKOV’s approach is analogous to a “‘determinate’’ approach 
in which the class of permissible initial data (see Section 1) is narrowed. We 
present here a statistical approach similar to the determinate one in which the 
notion of correctness according to TYKHONOV or of quasi-solutions is introduced. 
In this case we shall use the notions and definitions of information theory 
presented in the paper by DoBRUSHIN [49]. 

As it was noted in Section 1, most of the problems of mathematical physics 
reduce to the examination of equation 


Ag=f (1.1) 


where ¢, fare elements of the metric spaces ®, Fand A is a continuous operator. 
In practice the right-hand side is, as a rule, a result of observations made with 
the aid of real physical instruments and therefore cannot be considered to be 
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given with absolute accuracy. However, in practice it is often the case that a 
certain statistical distribution of errors in the definition of the right-hand side 
of (1.1), and the statistical distribution of solutions @ are given. Motivated by 
this we propose the following approach. 

Let ®, F be measurable spaces (sec [49), i.e., in the spaces ®, F besides the 
metrics, there are also given some a-algebras S,. S, of subspaces of the spaces 
@®, F and the operator A is measurable. In addition let 2 be a probability 
space, i.e., a measurable space with o-algebra S, and with probability measure 
P { } given on this o-algebra. The elements ¢, f in (1.1) will be considered to 
be random variables with their values in the spaces ®, F. Morcover, we consider 
also a random variable f with values in F—the right-hand side of (1.1) with 
an error—and the random variable (f, f) with values in F x F, the product 
of F with itself. Let B be an operator with domain of definition in F and the 
range of values in ®, and let @ = Bf, where gy, @ are random variables in the 
space @ x @, 

The problem of constructing an approximate solution to (1.1) consists of 
constructing such an operator B that the distance between gy, @ be sufficiently 
small where the notion of smallness should take into account both the metric 
of the space ® and the distribution of random variables (9, §). 

As an estimate of the deviation of the approximate solution @ from the 
exact solution @ one may take, for instance, the quantity mo(@,@)* or the 
quantity P {o (y, @) > eh. Corresponding to these, we introduce two notions 
of optimal solution of equation (1.1). 


1. The function @ = Bf is called an optimal solution of (1.1) if 
mo(g,@)=min. 
2. The function 6 = Bf is called an Sotiinal € solution of (1.1) if 
P {o(9, De eel 


If the above estimates of deviation of  from¢ do not exceed the number 6 
we call @ a solution of (1.1) with accuracy 6 in the sense 1 or 2, respectively. 
This statement is in its content close to the decoding problem of information 
theory (sec [49]). 

We note that since in our statement the operator A is determined, the 
distribution of the random variable gy in © uniquely induces the distribution 
of the random variable f in F while the joint distribution of f, f inFxF 
induces the joint distribution of the random variables 9, f, in the space ® x F. 
Since f in the statement is considered to be given, we may consider the dis- 


* mis the mathematical expectation, the random variable 0 (9, @) is assumed to be 
measurable. 
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tribution in ® induced by a joint distribution of g, f at fixed f. If the dis- 
tribution is to be given, the above formulated problem of constructing optimal 
solutions and solutions with prescribed accuracy become classical problems 
of the theory of probability. However, the practical construction of the joint 
distribution of g, es may involve considerable difficulties, and so in most cases 
it is expedient to choose simpler although less accurate methods. 

The determinate approaches in which the notion of correctness according 
to TyKHonov or of a quasisolution are introduced, may be regarded as parti- 
cular cases of the proposed approach. In fact, let the set M introduced in the 
previous section be measurable, and 


P{peM}=1. (1.10) 


Clearly the equality (1.10) may be regarded as a priori information that 
the solution @ exists and belongs to the set M. 

Proceeding from the above ideas, we now give some results relating to 
optimal solutions and to solutions with a prescribed degree of accuracy, 
obtained under some additional assumptions. 

1. Let F, ® be Hi.perr spaces. We denote by p za random variable in ® 
whose distribution is induced by a joint distribution of , f for fixed re One 
can easily see that in this case the optimal solution in the sense of 1 is given 
by the formula, 


6=JSp;dw (1.11) 
Q 
In fact, for any @ 


m|1G—pyll?=I16||?-2(G, Jpzde)+t fllysll? do 
2 2Q 


from which follows the validity of our assertion. 


2. Let the solution of (1.1) be unique, and let (1.10) hold, where M is some 
compact set and M, is the image (in the space F) of the set M. Then, as noted 
in the previous section, the operator A~! on a set M, is uniformly continuous. 
Its modulus of continuity will be denoted by « (1). 

Next let 


P{o(f,f)>5}=0 (1.12) 


where 6 is a constant. 
It can be easily seen that if ¢ in definition 2 of an optimal solution satisfies 
the inequality 
ex<a(d) (1.13) 


then the quasisolution of (1.1) defined previously with the set M is an optimal 
solution in the sense of 2. 

In fact, let @ be a quasisolution of (1.1) with the righthand side 4 and the 
set M. 
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Then, by virtue of (1.10), (1.12), (1.13) and the definitions of quasisolution, 
and modulus of continuity « (t), we get 


P {0(9,)>e}<P {o(f, f)>5}=0 
which means that the quasisolution @ is ¢ optimal in the sense of 2. 


3. Let B, be a regularizor of (1.1). 
We estimate the accuracy, in the sense of 2, of the solution of (1.1) with the 
aid of the regularizor B,. We assume, in addition, that the operator B, is uni- 


formly continuous on F for any 2 > 0, and denote by « (e, 2), o, (t), B (t) the 
functions 


a(e,2)=P {0(B,4Q, p)>e} 


a%,(t)= max Q(B,u, B,v) 
e(u,v)<t 


B()=P {oF f)>t} 


By virtue of the above assumptions and the known properties of the random 
variables 


P {o(BaF.)<0}=max| PY 0(B,7,B,0<St: Plo, ost | (1.14) 


Plows a,fr<3}2P}o.s)<e(5)}=1-a|a:(5)| (1.15) 


where «; (t) is the function inverse to «, (t). 
Substituting (1.15) in (1.13) we get 


P {0(B,f, 9) >e}<max EE (=)| 5 1)| (1.16) 


We note that the estimate (1.16) makes possible an optimal choice of the 
parameter A of the regularizor B,. 


4. We now give the estimates of the effectiveness of a particular regularizor 
for the Caucny problem for the LAPLAcE equation considered in the preceding 
section. Let F and @ be the space L, of functions defined on the segment [0, 77] 
and A, be an integral operator 


Anp= 1 Ky (x, £)p (2) dé 


with the kernel K,, equal to 
2 2 ‘ , 
K, (x, 3 me) cosh" kh-sinkx:sinké. 
0 


Next, let the distribution of the random variables g, f, f be such that 
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“= 


P{\|A,, ‘oll >t}=(p2t7 +1) * 


= (1.17) 
P{If-fl>ga@r+h 
Asa regularizor B, we consider an integral operator with the kernel 
¢ 2 — KR des : y 
0,(x, = ¥. coshkh(i+de")* sinkx-sinkg hy >hy 
TO 
Onc can show that the operator B, satisfies the following relations 
B,|| <1" -r, (hh 
[Ball 1 (h, hy) (1.18) 


An, pl] =A" «ry (hy, ha) (E— BA) ll 
where 
ry (h, hg) = hh” (hy ~ hyOr-P- hy 
r, (hy, hh.) =(h2- aN . hoe hy 
From (1.17) and (1.18) we get 


P{|B,F—/>e}sP {IB If —fil>eh=(Qarn” ery 7 41)¢,” 
P{||B,A—E) gl|>e} <P {I A, ol >a rye} = (Pn re, 
from which it follows that 


P {\|B,f —oll>2e}<max[(q7a7 e777 4:1) 15 (p22 7 1) ]. 
(1.19) 


It is evident that the first expression in the right-hand side of (1.19) mono- 
tonically decreases as A increases, while the second one monotonically increases. 
This means that in the sense of the estimate (1.19), the value of the parameter 
of the regularizor will be optimal when the two expressions coincide. By 
equating the expressions in the right-hand side of (1.19) we obtain 


ho/[(hath 
et DP oirs 2/(hthr) 
ce) q & 


P{||B,f —¢ll S Deh (pyr a gilt Pin gn eg. ed) ; (1.20) 


From (1.20) it follows that for any 6, e > 0 there can be found a sufficiently 
large g— the index of the ‘“‘quality of the measurements” —so that the estimated 
probability will be less than 6. 

We also note that the regularizor B, used by us also permits optimality 
with respect to the parameter /,, but the analytical determination of the optimal 
hy leads to extremely complex calculations. It seems that the parameter should 
be chosen on the interval [H, 2 #7]. 


Chapter II 


Analytic Continuation 


The problem of analytic continuation is investigated in this chapter. Esti- 
mates of “‘stability” for the solution of some problems are given, and concrete 
solutions are constructed. 


§ 1. Analytic Continuation of a Function of One Complex Variable from 
a Part of the Boundary of the Region of Regularity 


The problem was first considered by CARLEMAN in [3]. The results were 
obtained in [5]. 


I. The First Problem. 


Let f (z) be an analytic function, regular and bounded in some bounded 
domain D so that 
If(|<M zeD (2.1) 


I’ is the boundary of D; J” isa part of and J” = J’ — I”, Let the values of 
f (2) on I” be known, and suppose it is required to determine f (z) in some 
part of D. We now prove a theorem characterizing the stability of the solution 
of the problem. 


Theorem 1. Let f (z) on the curve J” satisfy the inequality 


[f(z)ise (zel’) (2.2) 
Then the inequality 
Fs? Oe) (2.3) 


will also hold where w (z) is the harmonic measure of the curve J” with respect 
to the point z and the domain D. 
For the proof we consider the function 


p(z)=In|f (z)I. 


It is known that ¢ (z) is a harmonic function. From (2.1), (2.2) it follows that 
the function @ (z) satisfies the inequalities 
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y(z)<lIneé ze!’ 
y(z)<InM ze” 

from which it follows that 
y(z)<w(z)Ine+[1—w(z)]InM (2.5) 


The inequality (2.3) to be proved follows from (2.5). 

The problem we are considering can be reduced to an integral equation 
of the first kind and can be solved by one of the general methods. In fact, we 
denote by D’ the part of D in which f (z) is to be determined. We denote by 
P, Q the ends of the curve J” and consider some points P;, 0; (J = 1, 2) lying 
on I” sufficiently close to the points P, Q, respectively. We denote by ra j the 
part of J” lying between the points P;, Q;, and let the location of the points 
P;, Q; be such that I’; € Ij (see Figure 1). 


(2.4) 





We draw a curve J’; through P,, Q, which lies inside D so that the domain 
D’ lies inside the domain D{, bounded by the curves /%,, I}; D’€ Dj. It is evident 
that in order to determine the function f (z) for any ze D’ it is sufficient to 
determine f (z) on the curve J’, in view of the Caucuy formula at ze D’ 


1 AG 
Ti por 6—zZ 
and the f (z) values given on I’. 
We denote by /, (z) the limiting values on the curve J°5 of the Caucuy-like 
integral of the function f(z) on the are I” as z tends to I’} from inside of D; 
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and by /f; (z) we denote the values on the curve J’; of the Caucuy-like integral 
of the function f(z) on the are J” 


fil=zo lim je () Dat, zelr,. 
besry' 
zeD 





7a f= 55 (ES 


The limiting values of the Caucny-like integral are calculated by the familiar 
SOKHOTsKY formulas, so that for an investigation of this problem we may 
suppose that the functions f5, f, are given. 

We denote by 4 the operator transforming the arbitrary complex function 
@ (z) on I’ to the function w (z) given on Ij by the formula 


oO a6, 


W(z)==— zel. 


a 


It can be easily seen that the validity of the equation 


A@(x)=f(z)-fia(z)  zel2 


is necessary and sufficient for the value of the function @ (z) to coincide with 
the values on J’, of the unknown function f (z). Thus, the problem of deter- 
mining the function f(z) on J’, is equivalent to that of solving the linear 
operator equation of the first kind 


Ag=f—fr2. 
Then by the Caucuy formula for any ze I’, 
1 
f@=zq [Oat n@ 


from which it follows that the function f(z) on J’, is equal to 


f (z)=BW(z)+f,(), 
where W (z) is a function defined on J” 
IW(@I<M, 


(M is a constant depending on D), and B is the operator transforming the 
complex function W (z) given in I’ to the function BW (z) given on I’, by the 
formula 


prea 





BW (z)= 


It is clear that B is a linear completely continuous operator and therefore 
in the problem of solving the equation 
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Ag=f—-hf, 
the conditions under which one treats linear operator equations of the first 
kind are met. 
We now consider another method of concrete solution based on the integral 
formulas of CARLEMAN type. We call a CARLEMAN function for the domain D 
and the curve J” any function G (z, 2, 6) having the following properties 


1) GG,6,5)=-—+ Ol, C,6) 


where G (z, ¢, 6) is an analytic function of the variable ¢, regular and bounded 
in the domain D. 


2) the function G (z, £, 6) satisfies the inequality 
f IG CA ld <s. 
By 


We will construct the concrete solution for the problem with the aid of a 
CARLEMAN function. Let the values of f(z) on the curve J” be known with 
accuracy €, i.e., a function f, (z) is known on the curve J” such that 


lf(2J—-fS(zise, zer’. 
We denote by the function f; (z) 


1 ‘ 
IsD=z 5 f G(z,¢, 6) f,(0) al, 
Ti jp 
and estimate the difference f (z) — f; (z). 
According to Caucny’s formula 


=x | GLAS Oat. 


Ti > 
Therefore 


fO\-hO=zq J OLAS Oat 


tra J GLO O-F(0] de. (2.6) 


In virtue of (2.1) and the second property of CARLEMAN functions the first 
integral in the right-hand side of (2.6) satisfies the inequality 


IJ Gf (QHde[<6-M. (2.7) 


In virtue of (2.2) the second integral in the right-hand side of (2.6) satisfies 
the inequality 
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I CLF O-£.0] dl|<p(z,6)-e-y’, (2.8) 
where y’ is the length of the curve J”, 
H(z, 6)=max |G]. 
CeD 


Substituting (2.7), (2.8) in (2.6) we obtain 


If ©) fy) S 5 [5M + ne, 6)-0-7']. (2.9) 


Let the CARLEMAN function be known. We denote the root of the equation 


6 ey’ 
CZ) a 5) M (2.10) 
by 6 (z, £), and put 
6=6(z,&) 
in (2.9). Then the inequality (2.9) becomes 
1 
If (z)— fs(2)| SM “8 (2, 6) (2.11) 


It is evident that as e > 0, 6 (z, €) and therefore the right-hand side of (2.11) 
tends to zero. Hence, if the CARLEMAN function G (z, &, 6) is known, the func- 
tion f;, (z) may be considered to be an approximate solution of the problem. 
The difference between the “true” solution f(z) and the approximate one 
fs, (z) is thus estimated by the above inequality (2.11). 

It is evident that the smaller the function y (z, 6) is, the higher will be the 
precision of the approximate solution obtained by the above method. We shall 
now display a family of CARLEMAN functions with minimal » (z, 6) in the case 
where D is a simply connected domain. 

Let 0, (z, €) be the function which is harmonic with respect to the variable ¢, 
regular in the domain D and on J’ takes the boundary values 


0, (z,6)=In|z—{I, Cer. 


We denote by 0, (z, ¢) the harmonic function conjugate with respect to 
0, (z, €) and by 6 (z, ¢) we denote the analytic function 


0(z,0C)=8, (z,E) + i02 (z, O —9; (z, z) — 10, (z, z). 


Let w (z) the harmonic function conjugate with respect to the harmonic 
measure @ (z), and W (z) the analytic function 


W (z)=@(z)+i@(z). 


2 Springer Tracts, Vol. 11, Lavrentiev 
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We consider now the function 


@;-(z, Ind 
G(e,6, =p Lemp Hee w(+0(2,0-a2,9} 


where 





q(z, jy Gee | MEDD a) +0, ea], 


It can be easily seen that this function is a CARLEMAN function, 


u(z, symexp| BG D+S of, (2.12) 


The constructed CARLEMAN function is best possible from the standpoint 
of estimating the precision of the approximate solution by (2.11). However, 
it turns out that for certain domains the construction of the function is some- 
what cumbersome. Therefore, in order to solve the problem it is preferable, 
in most cases, to use another CARLEMAN function, namely ,the one pointed 
out by CARLEMAN himself (see. ref. [3], [29]). This function does not give high 
accuracy, but its actual construction is rather simple. 


§ 2. The Cauchy Problem for the Laplace Equation 


The Caucny problem for the LAPLACE equation in the plane is similar to 
that of analytic continuation considered in § 1. An analytic function is te be 
determined from its values on a curve on which Caucuy data are given. 

In fact, let the value of the harmonic function u(x, y) and its normal 
derivative 0/On u (x, y) be known on some curve J*, We denote by f (z) (z = x 
+ iy) the function f(z) = u + iv where v is the function conjugate to wu. It is 
well known that on the curve I’ 


v(z)= J ou (z)ds+C 
where Z, is one of the endpoints of I’. 

Hence, if wu (z), 0/On u (z) are known on J’, one may consider that the 
values of the analytic function f (z) on I’ are known. 

This section deals with the extension of the CARLEMAN function method for 
the case of harmonic functions in a thre-dimensional space. 

Let D be some domain of three-dimensional space bounded by a sufficiently 
smooth and closed surface 2. Let 2” be a part of X, and 


2 Us =>. 
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We call a CARLEMAN function for the domain D and the surface Y any 
function of the three-dimensional vectors x, £ and the scalar parameter 6 
satisfying the following conditions 


I) G(x, ¢,6)=—— + G, (x, ¢, 6) 


eB 


where r (x, ) is the length of the vector x — é, and G; (x, é, 6) is a harmonic 
function of €, regular and bounded together with its gradient in the domain D. 


2) The function G (x, &, 5) satisfies the inequality 


| 
sa 516+ 


= 
S 








qt torso (2.13) 


where 0/0”; is the normal derivative with respect to the surfacc 2, and do; is 
the element of area on & with respect to the variable é. 

Let the CARLEMAN function G (x, &, 6) for the domain D and the surface 
be known. By making use of the function G (x, &, 6) we can prove a “‘stability” 
theorem and carry out the method of concrete solutions of the Caucuy problem 
for the LAPLACE equation in three-dimensional spacc. We dcnote by y# (x, 0) 
the function 


Theorem. Let the harmonic function u (x) be regular inside D and contin- 
uously differentiable over the closure of D, and let it satisfy the inequalities 


p(x, d)= max{ + | a =—G 

















luo u Os) <e xed’ 
= (2.14) 
Ju cor+}2we9 <M xed” 
Then the inequality 
ju (x)| <2 Mt (x, €) (2.15) 
also holds; where t (x, €) is a root of the equation 
u(x, t)-e=Mt (2.16) 


It is evident that t (x, €) 7 Oas e> 0. 

It is clear that the value of the harmonic function uw (x) at a point x inside 
D is expressible in terms of the u, 0/Onu values on the boundary D accor- 
ding to GREEN’s formula 


1 0 fe) 1 re) 0 
un=zz J \ogu- ane cou dept ae f {c aa Oa uf de. 


20 Analytic Continuation 


from which, by virtue of (2.13) and (2.14) it follows that 
ju(x)| <u(x,6):-e+M-d. (2.17) 


Assuming that 6 in (2.17) is equal to a root of equation (2.16) we obtain 
inequality (2.15). 

The construction of the concrete solution with the aid of the CARLEMAN 
function is completely similar to that given in Section 1. Let the Caucuy data 
for the function w (x) on the surface X” be known with accuracy eg, i.e., the 
functions @, (x), p, (x) are known on the surface 2” so that 


\p,(x)—u(x)| Se, 
xed’ (2.18) 


ve) Zu (x)] <e, 








and let the function (x) on the surface 2” be known to be bounded together 
with its normal derivative so that 


Ju (x)| + <M, xed”. (2.19) 








6) 
On u(x) 
We consider the function 
1 a) 
wO=zz J E v.O-3-G-9, o| dog, 


and estimate the difference u(x) — u; (x). Substituting the corresponding 
integral expressions for u, us and utilizing inequalities (2.13), (2.18), and (2.19), 
we have 


1 a) ) 
\u(x)—us(x)|= AS E ant! -Z6u] do, 


1 0 0 
Fae \¢| 2u-v.]-2< [u-ea} do, 


t 


S 
<2M6-+y(x,5)8-7— 


(2.20) 





where SS’ is the area of the surface 2”. 
Let 6 be a root of equation (2.16). Then inequality (2.20) takes the form 


4n 


Hence, if the CARLEMAN function G is known, the function u; (x) will be an 
approximate solution with accuracy estimated by the inequality (2.21). 


|u (x) —us(x)|<Mrt(x, o(2+2). (2.21) 
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We now tum to the problem of constructing the CARLEMAN function. The 
existence of the CARLEMAN function for an arbitrary simply connected domain 
D and for any part of the boundary 2’ bounded by a smooth curve follows 
immediately from results of MERGELYAN, [47]. In that work it was established 
that for any arbitrary pair of continuous functions @ (x), wp (x) defined on a 
smooth piece of the surface £’ and for any 6 > O there exists a harmonic 
polynomial H,, (x) of order m satisfying the inequalitics 


|H,,(x)—@(x)| <6 


<5. (2.22) 





Ha ()— vO) 





In this paper an algorithm is given for the construction of the poly- 
nomial. 

We consider the function [r (x, £)}-1, and let the point x lie inside thc 
domain D. We fix 6 and construct a harmonic polynomial P,,, (x, €) of order #7 
with respect to the variables ¢, satisfying the inequalities 


lr Gs =k (x, €)| <4 
fer” 


0 col E, fe) 
an! (x, 5) — a Pm $) <6. 


We denote by G (x, &, 6) the function 


G(x, €,5) =r +(x, €)—P,, (x, €). (2.23) 


It is evident that the function G defined by equality (2.23) is the CARLEMAN 
function for domain D and surface 2”. 

The above method for the construction of CARLEMAN functions is distin- 
guished by its great generality. However, the algorithm for such a construction 
is very complicated. We now give a simpler method for the construction of a 
CARLEMAN function for certain values of x. 

We denote by 2 (x, 7) the sphere with center at the point 7 and intersecting 
the point x. Let the point x € D satisfy the condition: There exists a point 7 
such that the surface 2” lies inside the sphere 2 (x, y). By expanding the 
function r—! (x, €) in TAYLOR series in the variables €, — 4,, we obtain 


PIG HS) wnt) (2.24) 


where P, is a homogeneous harmonic polynomial of order k in the variables 
& — n; (J = 1, 2, 3). We denote by G,, the function 


G, (x, 9, 6)=r7* (x, €)- 2 Pam e—n). 


22 Analytic Continuation 


It may be casily seen that the function G, satisfies the inequality 
20 \c—H\ 
Ix—nl"** (Ix—nl-1€—al)? 


In fact, on any ray € — 1) = 2 ty where 7 is a scalar and y is a unit vector, the 
sum 


|G,|+|grad G,| < (2.25) 


oy P(x, €-—m= a2 P,(x,,T, Y) 
0 0 


is a partial sum of the TAYLor series in the variable t for the analytic function 
1 


F(a,x4,4.))=—_—_—_———___—_—.. 
Sa ae Ter ea Fe 





The inequality (2.25) then follows from well know estimates for the reminder 
term for elementary analytic functions. 
From (2.25) it follows that if the number n satisfies the inequality 


9 n df 
CeCe ee eee (2.26) 
Ix—yl"™* [[x—n|-RO)] 
where S”’ is the area of the surface 2” and 
R(Z")=max|E—n], ed”, (2.27) 


then the function CG, is the CARLEMAN function for the domain D and the 
surface 2”, 


§ 3. Determination of an Analytic Function from its Values on a Set Inside 
the Domain of Regularity 


Let f(z) be an analytic function, regular in the unit disc D, 
If(z)|<1, zeD (2.28) 
and let A be some set inside the disc 
lz|<R<1 (Dp). 


The problem is to determine f (z) in the entire disc Dp from its values on A. 
It is evident that for the solution of our problem it is sufficient to determine 
J (2) on the circumference 


|z|=R (Fr) 


since according to the Caucuy theorem 
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PO=s5 f ah OMe. 2.29) 


We denote by L the operator which transforms the function given on I", 
to the function given on the set A by the formula (2.29), and by L, the operator 
which transforms, by means of the formula 


$O=sa Sal OM, 


the function given on J’; to the function given on J’p. It is evident that L, L, 
are linear, completely continuous operators. 
Solving this problem is equivalent to solving 


Lo=f (2.30) 


where f is the known function defined on A, andg is the function to be found 
on J',. and 


gp=Liy (2.31) 


where yp is a function defined on J’, satisfying the inequality 
lpl<t. (2.32) 


Hence, the solution of the problem is reduced to the solution of an integral 
equation of the first kind (2.30), i.e., to the determination of the function g 
from the known function f, provided it is known that the function ¢ satisfies 
the conditions (2.31), (2.32). 

Now, we prove a theorem which characterizes the stability of the solution 
of the problem. 


We introduce first some definition sand notation. Let A be aset of points on 
some curve. Let us denote by H, (A) a system of n intervals containing the 
set A 

H,(A)>A 


and by « (H,,) the linear measure of H,,. We call the infimum of 4 (H,) for all 
possible H,, containing A the n-measure of the set A and denote the m-measure 
of A by 

yi,(A) =inf (H(A). 


It should be noted that this m-measure is unusual since, in general, it does 
not possess the property of additivity. However, it is evident that m-measure 
possesses the following property which replaces in some sense the additivity 


bay (Ay) +4, (Ax) > fn, +n (A, + A) : 
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Now, let A be a set of points of the complex plane lying in the disc Dp. 
We denote by A, the radial projection of the set A onto the circumference 
|zj=r. 


Let W(z) be a conformal mapping carrying the unit disc onto itself, and 
carrying the set A onto the set A,, < D,, where D,, is the annulus 


0<r<|z|<o<1 (D,.) 
By yf? (A) we denote the supremum of x, (A,,,) taken over all possible 
mappings W (z) 
Hy (A)=sup fy (Ayr) - 
Theorem. Let A be a set lying inside the disc 
lzZ|\<R<1 (Dr) 
and let the analytic function f (z), regular in D, satisfy the inequalities 
lf (z)| <1 (zeD) 
If (Zl Se (ze A) 


Then the inequality 


Ine(_l—|z 
\f (Z| <exp \- Cy ine—leD} (2.33) 
In p72 (A) 
is valid, where C, is some constant dependent of the radii 9, r, R and p satisfies 
the inequalities 
weAP [wea Po 
a ey eae 
| C5 sexs C; (2.34) 


with C, also dependent on g,r, R. 
First we prove two lemmas. 


Lemma 1. For any points z, a lying in the annulus D,, the following ine- 
quality is valid 


a-Z a,—Z 


1—4,z, 


r 





1—az 














where a,, z, are the radial projections of the points a, z on the circle I’, 
|jz|=r. 


Without loss of generality we may consider the number a to be real and positive. 
Let first |z| > a. We consider the function 


a—TZ 2 





I,@= 








1l—atz 
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Differentiating the function /, (t) we obtain 


br qyag tC? +97) +a? ax) 
, |1—a-z|* 





ms yee 
Let 


x<a 
then 


(x? +y?)(L+a?—ax)—ax>x?4+y?—a?>0. 
Now let 


x>a 
then 


(x? +)?)(L+a*—ax)—ax>(x—a)(1—ax)x>0. 
Hence 
I, ()>0 


from which it follows that for |z| > a, 


a-Z aA—Zjq4 








oe (2.35) 








1=4az;,; : 
Now, let 
|z| =a]. 
We consider the function 
2 


L@=t|\-= 











1—taz 


Differentiating the function J, (t) we obtain 


p74 ea) ty Ja") 
1,()= (—a?)* >0 





from which it follows that for |z| = |a| 


a—-Z 


a,—Z, 
1—az 


(2.36) 














1—a,z, 
Finally, for |z]|< |a], we consider the function 


2 
T—Z 


I,(t)= 











1—tz 
Differentiating the function Jz (t) we have 


_j[a?Gtty)=2ax+ 1] (a—x)- (0? +y?-2ax+a7)(a-[2l?—¥) 


e jl —az|* 
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Let 
a(x*+y7)—x<0 
then, clearly 


I,(a)>0. 
Now, Iet 
a(x*+y7)—x>0 
then 
a—x>a(x’+y’)—x 
and, hence 


Ba 2 ee: 
—-xjqad- 1- 
peje eG Se 

\l1—az*| 
from which it follows that for |a] > |z| 


a-—Zz Qj) —Z 





(2.37) 














l1—az 1=a,-\*2Z 


The inequality to be proved follows from inequalities (2.35), (2.36), and 
(2.37). 
Lemma 2. Let 4,, a2, Z,, Z, be some complex numbers, and 
la,|=|a2|=|z;]=|z.;=r<1. 
Then, if the inequality 


a a 
arg—<arg—<n 


is true, then the inequality 
| 


Z4—ay, Z2—@s 


1-—@,z, 














1 | a421 
will also be true. In proving the Lemma, one may suppose that 
a = az =a 


where a is real and positive. 
Let 


Z;Sae"": .2z;=ae"*. 
We consider the function 
24a |* 
l1-—az 





I(g)= z=ae 








Differentiating the function J (g) we obtain 
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28 2\2 
1, . a sing (l—a’*) 
I amy | ea >0 


from which we obtain the inequality to be proved. 


We now prove the theorem. Let 1 = w (z) be the conformal mapping of D 
onto itself for which 


Hp (A) = fy (Ayr) 
Under this mapping, the set A is carried into the set 
A, Cc Dy: 
Let us consider the function 
y(z)=In|f (z)|. (2.38) 


It is known that @ (z) is a harmonic function, regular at those points of D 
where the function f (z) is not zero. We consider some point a belonging to A,,. 
According to the condition of the theorem 


If (a)|=e,< (2.39) 


From (2.38), (2.39) and well known properties of harmonic functions it follows 
that 


If (Zl =e, 


on some analytic curve y, intersecting the point a (if f (a) = 0 then the curve 
consists of the single point a). 
Let the curve y, be an open arc in the annulus 


2 1 2 
—r<|zl<— 2 D 
zr Slel<z (1420) 


te 


Then the curve y, either intersect one of the two circles 
2 z 
|z| =a F Ci r) 
3 


1 
lel=z+20)  (f,) 
or bounds a domain containing D273 ,. 


Suppose 7, intersects the circle F , in some point a,. The distance between 
the points a, a, satisfies the inequality 


la-alear. (2.40) 
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It is not difficult to show that by virtue of (2.4) the harmonic measure of the 
curve y, in the domain D with respect to the point z satisfies the inequality 





r 1-|2| 
SS 
O12, We TH Iz] 
from which it follows that 
r 1-|2| 
<—_ —__] 2.4] 
OOS ra Teal ey 
From (2.41) we obtain 
r 1-|2| 
< ] 2.42) 
el exp {nee at 


In case the curve y, intersects the ring 
1 
|z|==-(1 +20) 


we similarly obtain the inequality 








1-—@ 1-|z| 
Lf @|sexp fine mn iolelt (2.43) 


Finally, if the curve y, bounds a domain containing D;3, then from the 
maximum modulus principle, in the whole domain D 2/3, the inequality 


If@l<e  zeD2, 


is valid, from which the assertion of the theorem immediately follows. 

Thus, if just one of the curves y, is not closed in the annulus Dios t the 
statement of the theorem follows. Now suppose all the curves y, are closed in 
DD, and consider any curve y,. As y, is closed in D, it follows that inside the 
domain bounded by the curve y, there is at least one zero of the function f (z). 
We denote by y the sum of sets y,, and by y, the projection of y on the circle 


lear. 


In the ring 


ZrsiszQto) Bb 


te 


suppose the function f (z) has p zeros at the point a, (k = 1, ..., p) (roots are 
counted according to their multiplicity), and let 


F(z)=f (2): = — HE 
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We prove now that in D,, there exists a point z* such that 


1 Pp 
1a ——_—_— 
IF ( =| | (2.44) 


Let z be some point of D,.. and z, be the radial projection z onto Fe 
By virtue of Lemma 1! we have the inequality 


(2.45) 














Let us consider now the set y,. Since the function f (z) has p zeros in Diss 
it follows that the set y, consists of no more than p intervals. Then by virtue of 


Vr? Avy 
it follows that the measure of y, satisfies the inequality 
UO) = Hy (A). 


We denote byy/ the part of y, which belongs to the semi-circle 0 < argz< 2, 
and by y” the difference 


r=), 


It is evident that one of the two inequalities holds 
ee eee 
LODZ HO) 


itd 1 
i (yy) 25 i (y,) . 
For the sake of definiteness assume 
re | 
LOD2> HO). 
By virtue of Lemma 2, in proving (2.44) one may consider y, to consist of 


one interval. From CHEBYSHEFF’s results on polynomials which deviate least 
from zero it follows that on the interval y} there is a point z* such that 


[Tlaee— 241 [Caer > [ue (aye, (2.46) 


where C, is some constant which depends on r. From (2.46) and (2.45) we see 
that on the set y we can find a point z* such that 
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* 1 Pp 
(2.47) 
5 Cup al 


1 a,—z* 


Ls 1 —a,Z 














and the inequality to be proved follows from (2.44). 
Now we consider the function 


v(z)=In |F (z)| 


The function F(z) has no zeros in the annulus D, » and therefore the function 
vu (z) is a regular harmonic function in D,, . Then by virtue of (2.44) and the 
conditions of the theorem 


v(z*)< —pln pf (A)+Ine 


v(z)<O = (ze D,,) ve) 


Wc denote by cS the boundary of the annulus D5. It follows from (2.48) 
that 
J v@ds<C;, [Ine— pln py’ (A)] (2.49) 
Tig 
where C,, is a known constant dependent on the radii r, g and it follows from 
(2.49) that the function v (z), in the annulus D,,, satisfies the inequality 


v(z)<C,, [ne— pinp’ (A)] (2.50) 
From the inequality (2.50) we obtain 
C"re 
If (ls 5 re - + (2.51) 
Lup (A)]? 


It is not difficult to obtain from the inequality (2.51) the following estimate for 
the modulus | f (z)| in the whole disc D 


; C"re (1- I= 1) 
If (2s WONG (2.52) 


Then, since the function f (z) has p zeros in the annulus D ro it follows that 


If (Disv?, (2.53) 
where 
_lzl+6, . 2+¢ 
L+lzjo° =" 3 
It is easy to see that 
Inv<C,:(1—|z|) (2.54) 


where C, is a constant dependent on @. 
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Let p satisfy (2.34). Then by virtue of (2.52), (2.53), and (2.54) the inequality 


¢ xpd —C, (L—|z|)) 8 e 
F @isexp Ci Dee | 


is valid, and this concludes the proof. 


We cite some examples of obtaining ure (A) and the corresponding esti- 
mates. 


1. Let the set A,,, have positive JORDAN measure # (A,,,). Then, it is evident 
that 


iy (A) 2 a (A,,) ; 
Ky (A) rahe (A,,,) : 


nO 


If, in particular, A,,, Consists of 7 curves, then 


Uy (A)> Ut (Avr) n<m, 
in (A)=#(A,,), n>=m. 
If m1 = 1, then the estimate of the theorem, accurate up to the constant C,,. 


coincides with the estimates for an analytic function using harmonic measure. 
2. Let the set A consist of the points a, (kK = 1,..., 7) 


la,|=a. 


arga,=ko, (==). 


One may easily see that in this case 
i, (A,)=a(n—p)o, psn, 
U,(A,)=0, p>n. 
The relations (2.34) in this case are as follows 


[(n—p)d5-a]?<e°*<[(n—p+1)5-a]}? (2.55) 


It is not difficult to show that the number p satisfying (2.55) satisfies the ine- 
quality 
nd 
p = C In e’ ea Sap 
1+06'Ine 


from which we obtain 


no’ Ine’ 
<exp 4 C, ———_-(l-- -p} (2.56) 
If (z)| P| ine | 
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3. Let the set A consist of the points a, (kK = 1,..., 00) 


T 
la,|=a, CE 


One can see in this case that 
a 
H,(A)=—, (2.57) 


and the relations (2.34) are as follows 


GN Zeee an le 
2P ~ a p> 





from which 


p*>C7,/|Ine| 
and 
If (z)| <exp {—C,./ In] el-(1—[z))}. (2.58) 


§ 4. Analytic Continuation of a Function of Two Real Variables 


Let f (x, y) be a function of the two real variables, analytic and regular in 
the disc 


x*t+y*<1, (D) (2.59) 


and let A be some set inside D. Let the values of the function f (x, y) be known 
on the set A and suppose it is required to determine f (x, y) inside D. We 
prove a theorem characterizing the stability of the solution for a certain class 
of sets A. 

First some definitions are introduced. Let us divide the x, y-plane into a 
system of parallel strips H, of width h, (k = +1,..., + 00). Let T denote this 
division while A, denotes the part of the set A that belongs to H, 


A,=A0 H, 


By A, we denote the projection of the set A, onto the center line of the strip 
H,. Let t> 0, C;(j = 1,...) be some constants and let n, be an integer 
satisfying the inequalities 


A ny + 1 yz Mk 
Es mat 2| <t+Ch,.< Eo (2.60) 


where 44, (A,) is the n-measure of the set A,. 
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We denote by H,, (A, T) a set of strips Hi, for which the following inequality 
is true 


In(t+C,h,) 
Cyn, AD” ; (2.61) 


and by 44,2, (A, T) the m-measure of the projection of the set H,, (A, T) onto 
some straight line perpendicular to the lines of the division T. We call the supre- 
mum of “,,,, (A, T), for all possible divisions T, the (m, 2, t)-measure of the 
set A and we denote it by 


Lonar (A) = amp ee (A, T)} : 


Theorem. Let A be a set of points of the x, y-plane lying in the disc 
x?+y7<R*<1 (Dr) (2.62) 


and let f (z;, 2.) be an analytic function of two complex variables z,, z, regular 
in the hypercylinder 


(Rez,)*?+(Re aN a 
\Imz,}<%, |Imz,|<x. (2) 


If the function f (z,, z,) satisfies the inequalities 


If (21, 221<1, (z1,Z2)€2, 














(2.63) 
fy} <e (x, yeA, 
then the inequality 
If yy sexp 4 "4_| ox, yep (2.64) 
Sere Nada (Ale ee 
is valid, where the numbers m, / satisfy the relations 
Un+1,2 (A) dais -2 Linas (A) i 
LSE Ba Schaal < <= |-_—_—— 2.65 
C, Be Oa Gg. 2S) 
and C; are constants depending on R, x. 
We denote by 7, the division of the plane x, » for which 
Lime (A) —Einte (A, T)| <a. (2.66) 


One can obviously consider the corresponding strips H, to be parallel to 
the x-axis. Let the center line of the strip Hj, be the straight line 


Y= Vx 
It is not difficult to show that from (2.63) it follows that 


3 Springer Tracts, Vol. 11, Lavrentiev 
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[grad f(x, y)I<C3, (x, y)eDrp (2.67) 


where C3- depends on R and x. From (2.67) and (2.63) we find that at the points 
of the x, y-plane belonging to the set A, the function f (x, y) satisfies the ine- 
quality 


lf (x, yi<e+5 C3h,. (2.68) 


Let us consider the function 
~y,(zZ)=Sf (2; Vx), (x, y)EH,. 


From the condition of the theorem and from (2.68) it follows that the function 
Y, (z) is an analytic function, regular in the domain 
[Im z| <x; 
and satisfies the inequalities 
lp, (z)|<1, zEQ;,; oes) 
lp, (z)|<e+Cs3h,, (Rez, y,)€A,. 


Applying the results of the preceding section to the function @ , (z) we find that 
the following inequality is valid 


IRe2|<- 4" 





In(e+C3h,)-C, 
; < SS 
I, (z)| Scxp | In de (A,) > 
x 
Im z|<— 
\Im z1< 3 
where the constants C3 and C, depend on R, x hence 


In(e+C3h,)°Cy 
In pt, (Ay) 

Now we consider the function 

Wx(z) =f (x, 2). 


From the conditions of the theorem, the definition of the set H,,(A, 7) and 
(2.68) it follows that the function w, (z) is analytic, regular in the domain 


[Re z|</f/1—x? 





I Conisex | i (x, y)EH,. A Dp (2.70) 


(Q)) 
|Im z| <x 
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and satisfies the following inequalities 
lW.(zI<1, zEQ", 
(2.7!) 
lW.(yi<e?, (x, y)EH,,(A, T,) 9 Dr. 
Applying, as previously, the results of the preceding section we obtain, by 
virtue of (2.71) 
(xv, Rez)ED,p 


Ee re nee ee 
Yr. (2)I Sexp tr Hin2e (A, ot 


x 
Im z|<— 


from which it follows that 


2 . Cs 
Xx <exp | ———-——7.. x, yeDdD 2.72 
FG is p fate at (x y) R ( ) 
Passing to the limit at « > 0 in the inequality (2.72) we obtain the inequality 
(2.64) to be proved. 

Now we give two examples of the applicability of the theorem. 


1) Let A be a closed set with positive plane measure st (A). We show that 
in this case when 





he JInz|-C 
‘Sine (ADI 
the inequality 
1 
Linde (A)>oH (A) (2.73) 
is valid, and consequently 
Ine 
< - 2.74 
I vsexp| ont (2.74) 


where C is some constant. 
Let N be some number. We perform a division T of the plane into strips 


Ksysa (Hp 
and denote by H,,; the square 
bl ek 
N ~~" N 
(Hi,;) 
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We denote by p, the number of squares H,, of the strip H, containing 
points of A, and by p the total number of squares containing points of A 


lea} 
p=> Dy. 


One can easily see that the inequalities 


p>(A):N? 
(2.75) 

a Pe 2m 
Hn ( y= 
are valid. In fact, the first of the inequalities (2.75) follows from the definition 
of plane measure of the closed sets. In order to obtain the second of the ine- 
qualities (2.75) we note that by virtue of the definitions of p, and A,, on the 
Straight line 





there are p, non-overlapping intervals of length 1/N containing points of A,. 
Thus the sum of the lengths of any intervals covering A, cannot be less than 
the number 
Py-2m 
N 


from which the second inequality follows. 
Let us consider the strips H, for which the inequality 


1 
Pr2—q N-n(A) (2.76) 


is valid. 
Let (2.76) hold for 


k=k, (j=1,...,q) 


It is not difficult to show that the number q satisfies the inequality 
1 
q2q N-u(A). (2.77) 
In fact, assuming the opposite we see that 
1 2 
p<2Nq+(2 N-q) 7 N-n(A)< N*+ (A) 


which contradicts (2.75). 
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We denote by H the set of center lines of the strips H,; (j= 1, ....q), and 
by H the projection of H into the straight line x = 0. 

In analogy with the second of the inequalities (2.75) we find that for any a, 
the -measure of H satisfies the inequality 


q—2n 
N 





u,(A)> 
from which, by virtue of (2.75) and the definition of j,,,, (A, 7), we obtain 


n<ag N-a(A) 


Haac(A, TE HCA) (2.78) 


4S Int-C 
“In w(A) 


Passing to the limit as N— © in (2.78) we obtain the inequality (2.64). 


2) Let the set A consist of points p;; with coordinates 


- 2 i= | 

Na Sp T=1,...,@ 
4 ee: 

GTi j=i,it+l,...,.00 


We show that in this case for A < Cy Vint the following inequality is valid 


1 1 
Umar (A)= 5, 2V In} ( ) 
and consequently, 
If (x, yi sexp {—C |Ine|"/*} (2.80) 


Let us perform the division T of the plane x, y into the strips H, 


k-1 k 
2” on 
where N is some number. 
It is easy to see that for k = 2? (p = 1,..., N— 1) the m-measure of the 
set 4, is equal to 
= 1 
Ly (A=— (2,81) 


an tp es 


Therefore the relation (2.60) in this case are as follows 
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n+l Nk 
a <n <( c) (2.82) 
gmt ptt git P 








where 


From (2.82) it follows that for 


p< J|inz,|, 4<C,,/|In7,| 


the straight fine 


re a 
YrYer aN 





: (k =2°) 


belongs to the set A,, (A, 7). 
By virtue of the above statement, one can show that the measure of 
the projection of H,, (A, T) onto the straight line x = 0 satisfies the inequality 


1 1 
tm L43,.(4, T) | =— - —— 
Han Lb ( 1 Qm 2 [In tif 
from which follows 
1 1 
} ae | ) am QV Jin za ( ) 


Passing to the limit as N — 00 in (2.83) we obtain the inequality (2.79). 


§ 5. Analytic Continuation of Harmonic Functions from a Circle 


Let u (z) be a harmonic function regular and bounded in the disc 
|z|<1 (D) 
and let its value be known on the circle 


jzj=r<1 (C). 


u(re'?)= f (9) 


where f (~) is a known function. It is required to determine the function u (z) 
inside the disc |z| < r from the values of the function on I’,. Let us prove a 
theorem characterizing the stability of the solution. 
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Theorem 1. Let the considered function u (z) satisfy the inequalities 


1 . 
a J u?(re”) dose, 


(2.84) 
| 2 is ae 
— “(e*)do<1. 
3 a u(e")do< 
Then the following inequality is valid 
oe ae Ine-Ing 
pe + do< Se 
- ju (oe'?) ¢ ose lap } (2.85) 
It is well known that 
u(oe)= ¥ eX (a,coskp+b,sin kp)+ 2. 
1 
Let us write 
+ a 
Jagt be =C,. [Bae 
From (2.84) we have 
PC <e, 
(2.86) 


Let us estimate the function 
1 
‘i 


w= ¥ Ci =— J u? (ee) do 
0 Tt 


using (2.86). 
By virtue of (2.86) the sum to be estimated is bounded. Let the sum reach 


its conditional maximum over the variables C, for C, = C, (k = 0, ..., 00). 


It is evident that 
C,=0, k+p,@, 


where p, 4, Pp < g, are some numbers, since otherwise it is quite clear that for 
C,, one may choose variations 6, such that 


» (€,+6,)°e"*= > Cra: 
0 0 


D 


yi (Cp +5) < ¥, C2, 
0 0 
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Therefore C 3 C, satisfy one of the three relations 


C2+C7=1, ho 


Cor? PEC rE: 





(2.87) 
C,=0; (11) 
c: =(. (III) 
Let the first of the relations (2.87) hold. Then 
2q 
FA2 = —E—-r d 
Co pe??? 
2.88 
ze 2 p2P_¢ ( ) 
pee r74 


We note that from (2.88) follows 
rt<e<r??, (2.89) 
By virtue of (2.87), (2.88) we obtain 


re: pe Oe tee Pee (2.90) 
cs eR p2P__ p24 peep? : 


It is easy to show that for p < q the functions 
2p_ 24 


Q 
Fi (2. D=3 


2q° 
t P__4.°4 





2p,2q_ 4.24 ,2p 
@ —r'e 


F, (p, qQj= 


2 
per p24 


increase with increasing p as well as with decreasing g from which, by virtue 
of (2.89) and (2.90), it follows that 


oo sks : : Ine-In 
2 Cre"s lim _AeF1 (P. 9) + Fa (p, 9} =exp jnetneh (2.91) 
p.qr—s- 
Inr 


Now let the second of the relations (2.87) hold. Then 
3 C29? C2924 
. q 
0 
and by virtuc of (2.86) 
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from which we obtain 
ae ; Ine:l 
y C207 <exp {ineinel (2.92) 
0 


In case the third of the relations (2.87) holds, incquality (2.92) can be ob- 
tained in a similar manner. 

The statement of the theorem follows from (2.91), (2.92) and from the 
definition of the numbers C,. 

With respect to the function u (z) considered in theorem 1, let it be known 
that at the boundary of the domain the n-th derivative is bounded in ZL, i.e., 


a awe “a <M 
i eee dQ" Ps! 


and suppose it is required to determine the function u (z) on the circle 0 = | 
from the values of the function #(z) on the circle g = r. We shall prove a 
theorem characterizing the stability of the solution. 


Theorem 2. If the function w (z) satisfies the relation 


j u? (re’*) do<e; 





n 2 
| a u «| dg<l; (2.93) 
™ 0g” 


rt 


J u(e’)do=0, 


Tt 


then the following inequality is valid 
ee ee 1 
= J u°(e*)do<—, (2.94) 


where 7 is a root of the equation 


2Alnr—2nInA=lIne. (2.95) 


We note that for a sufficiently small ¢ the root of equation (2.95) satisfies the 
inequality 
2 1 Ine (2.96) 
4 Inr’ e 
The proof of theorem 2 is similar to that of theorem 1. From (2.93) it follows 
that 
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3 ye? Se: 
2 (2.97) 
y ReCi S15 
0 


where C,, as in the previous theorem, are expressed by the expansion coeffi- 
cients of the function u (ge) in a Fourier series in the polar angle g. 
In analogy with (2.87) we find that the numbers C,, for which the sum 


Gees i u’(e'*) do 
0 Tl -x 
attains its conditional maximum, satisfy the relations 
C,=0, k+p,q, 
while C,,, C, satisfy one of the three relations 
a Oe ae Oe C,=0 (ID 


pcr +q"Ci=1 c=0 (il) 


q 


from which the theorem follows in an analogous manner. 


§ 6. Analytic Continuation of Harmonic Functions with 
Cylindrical Symmetry 


Let x, y be vectors with components (%,...,X,), O44) ---> Yt) and let 
(x, y) be a harmonic function regular and continuous in the domain 


|x| <1, lyl<1 (D) (2.98) 


Let also the function u (x, y) as a function y depend on the variable r = |y| 
alone. It is required to determine the values of the function u (x, y) inside D 
from its values on the set P 


x yer, jl p=0: |x]|<o<1. 


We shall prove a theorem characterizing the stability of the solution. 


Theorem. If the function u (x, y) to be considered satisfies the inequalities 


lu(x,y)]<1, 9 (x y)eD; 
(2.99) 
lu(x,y)l<e, (x, y)eP 
then the following inequality also holds 
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lu(x, y)|<Ce’, (x, y)eD,, 


where D, is any domain lying strictly inside D and C, y are some constants 
dependent on D, and o. 

By virtue of the conditions of the theorem, the function « (x, y) may be 
represented as a function u(x,r) where r = |y|. The function (x, r) will 
satisfy equation 


0? 0 
ait = +4,u=0. (2.100) 


From the known properties of harmonic functions it follows that « (x, r) 
is analytic with respect to the variables x, r evcrywhere inside D. 

We will estimate the modulus of the function A?« (x, 0). We consider the 
function f (z) where z is a vector with complex components z, = x, + i&; 
(A = 1, ...,) analytic with respect to the variables z, and coinciding on a real 
hyperplane with the function w (x, 0) 


Ff (x)=u(x,0). 


[It is evident that the function f (z) is regular and bounded in the domain R 
of the -dimensional complex space 


zeER, if |Rez|<oe, |Imz|<h, o<e,<1, 
where / is a constant dependent on 0,. Define R, to be the domain 
zER, if |Rez|<e, |Imz|<h,<h 


From (2.99) as well as from the results of § 1 it follows that the function 
f(z) in the domain R, satisfies the inequality 


If (21<Cye"; (2.101) 


where C;, y, are some constants dependent on 9, /, @,, 4, whereas from (2.101) 
and from known properties of analytic functions it follows that 


If (1 sklate™ (2.102) 


where « is a constant. 
From (2.102) we obtain 


APu (x,0)<(2 p)!aje"! (2.103) 
We now expand the function wu (x, r) ina MACLAURIN series in the variable r 


Shes @ 1 o2P ( 0): .2p (2 104) 
u(xr)= ) Qpiaze x,0)+7 a: 


(the coefficients of the odd powers of r are zcro by virtue of the symmetry). 
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Substituting (2.104) in (2.100) we find that the functions 02?/dr2? u (x, 0) 
satisfy the following relations 
Dp aoe 


a ao (p=1,..., 0) 


07? 
eT Sera 


from which 


2p 





ieee Oe 


or2P say (x,0) (2.105) 


The validity of the statement of the theorem for any neighbourhood P 
follows from (2.104), (2.105) and (2.103). The validity of the assertion of the 
theorem for the arbitrary domain D, follows from the validity of this assertion 
for the neighbourhood P and from general estimates for analytic continuations. 


Chapter Lif 


Inverse Problems for Differential Equations 


We consider two types of inverse problems for linear differential equations. 

The first kind is the so-called inverse problem of potential theory where the 
equation is non-homogeneous, the coefficients of the differential operator are 
given and it is required to determine the right-hand side. In this case we shall 
restrict our consideration to the LAPLACE operator, although a number of 
results are carried over to more general elliptic equations. 

The second kind is that in which the equation is homogeneous and it is 
required to determine some of its coefficients. A well-known problem of this 
kind is the inverse problem for the STURM-LIOUVILLE equation. 


§ 1. The Inverse Problem for a Newtonian Potential 


By an inverse problem for a NEWTONIAN potential we mean the following. 
On a part of the surface of the unit sphere D, there is known a harmonic func- 
tion u (x) which is an external NEWTONIAN potential of a body with unit density 
which is star-shaped with respect to some point 7 and lies inside the sphere 

|x|<R<1 (Dp) 

It is required to determine the dimensions and shape of the body from the 
function u (x). This problem has an important application in geophysics and 
a number of papers are devoted to its investigation. 

The uniqueness of the solution of the problem was proved by Novikov. 
The notion of correctness in a compact set was introduced by TYKHONOV in 
an application to a problem. In [16] estimates were obtained which characterize 
as properly posed the problem for a class of domains bounded by continuously 
differentiable surfaces with bounded derivatives. 

Let the polar coordinates of the vector y — 7 be denoted by 0, @, @ and let 
the surface of the body inducing the potential u (x) be defined by the equation 


e=f (9,9). 


It is well-known that the function u (x) is equal to 
n/2 x f (¢,9) o’ sin fa) 


= =~ dodo dO 3.1 


where r (x, y) is the distance between the points x, y. 
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Equality (3.1) may be regarded as an integral equation for the unknown 
function f (¢, 0) which defines the body to be found. If the integral operator 
in the right-hand side of (3.1) is denoted by A then (3.1) will take the form 


u(x)=Af (9,0). (3.1’) 


The operator A in the right-hand side of (3.1’) clearly transforms a function 
of the polar angles g, 0 to a function defined on some part of the unit sphere. 
Thus, the solution of the inverse problem of potential theory is equivalent to 
the solution of the nonlinear integral equation (3.1) or (3.1’). 

Now we proceed to the presentation of a stability theorem for the inverse 
problem of potential theory. For simplicity we restrict ourselves to the case in 
which the part of the surface D on which the potential] « (x) is known, is defined 
by the inequalities. 

x,>H, |H|<1 (S,) 


Theorem 1. Let the functions fi (9, 9), f2(@, 0) satisfy the following con- 
ditions: 

1) The functions f(g, 9) (7 = 1, 2) are continuously differentiable for all 
values of the polar angles g, 9 and their derivatives satisfy the inequalities 





la @ 
Sie.o]<ms neo <M (n*—6). (3.2) 

2) The following inequality holds 
[uw (x)—u2 (x) Se, XE Sy (3.3) 


where 
uj (x)=Af;(Q, @). 
Then the difference f, (¢, 0) — f,(Q, 9) satisfies the inequality 
Cc 
If: (9, 8) — fa(, 0) <— (3.4) 


where Cj is a constant depending on M, and n is a number defined by the 


relations 
1 2nt+3 1/C3 1 2n+1 
(4) <C, <(+) (3.5) 


n 
while C; are constants depending on R, M. 
We first prove a lemma. 


Ine 











Lemma 1. Let w (y) be a summable function which does not exeed one in 
modulus and let the function f(¢, 0) satisfy the conditions of the problem as 
well as the inequalities (3.2). Then if the function 
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n/2 x f£(p,0) 2 


eo sind-@(y) 
u(x)= ———— doe dg dé 
©) a. J 0 r(x, y) eee 


satisfies the inequality 
ju(x)|<e (xeES,,) (3.6) 


then the following inequality also holds 
1 


ju (x)|+jgrad u(x)|<C, [Ine] <s (3.7) 


where the constant C; depends on R, M and x is any point lying outside or on 
the surface 


o=f (9,0). 
It may be shown from (3.6) that the function « (x) satisfies the inequality 
lu(x)j<e=e, ([x]=) (3.8) 


where Cg is a constant depending on the radii R and H. 
Let the polar coordinates of the vector x— 7 be denoted by 0’, y’, 9’ and 
the distance between the points x and y by 


a Ce Q’, Os 0,0, 0)=r(x, y). 


Let us consider the function 


Tn 2 BLO) 2 sin Ow 
CG an ae g 0) 


——_.——— dodo do. 
—nj/2—-n 10) r(Z,@ , 0’, 0, 9,0) ¢ 


It can be easily seen that w (z, 9’, 0’) is an analytic function of z which is 
regular in the domain 


Rez—C,|Imzi|<f(9',0') (Qy6) 


where C7 is a constant depending on the constant M in (3.2), and satisfies the 
following inequalities 


lo (z, Q’, 0’) | - |grad w(z, ¢’, 0”) | < Cg, ZEQ,. 9: > (3.9) 
d ’ 9 ‘ne 
Rez>f (qg’,0’); Imz=0; 
|o(z, @'; O’)| se, > 
(3.11) 


Rez>1, Imz=0, 


where C; are constants depending on M. 
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From the inequalities (3.9), (3.11), some properties of analytic functions 
and the results of Chapter II, it follows that the function « (z, 9’, 0’) satisfies 
the inequality 


lo (z, 9’, 0’) + grad w (z, 9’, 0')| <e°" 


12 
Rez>1, Imz<—C,. 4) 


By virtue of (3.9), (3.10), and (3.12) the functions w (z,q’, 0’) and grad w 
satisfy the conditions of Carleman’s theorem ((3] or [29]) and therefore the 
following inequality holds 


|o(z, 9’, 0’) + |grad w(z, 9’, 0')|< C1, exp {Ine |z—f (9', O17}. (3.13) 
Using inequalities (3.10) and (3.13) we obtain 


[u(x)| +lgrad u (x)| =10(0', ', 1+ 

+|grad o(0', 9’, O'S Cy, exp {Inela—f (g’,0')]7} + (3.14) 

+C6[o—f (9, 0')]- (1+ In Lo—-f (9, 0]; 

f(@',0)<o'<@<i+y. 

Setting 
1 
O—f (9,9) =|Ine| 2 C12 
and using the obvious inequality 
Leo 
[In [@—f (9, 8’) I< Cy3 Le—-f (9, 0] 


we get the inequality (3.7). 


Lemma 2. Let w (x) be a harmonic function regular inside the surface 
Ix—nl=f(g',e)+h, h>o (24) 


where f (@, 9) satisfies the conditions (3.2) and does not exceed one in modulus 
in this domain. Then on the surface 


Ix—nl=f (9', 6’) (2) 


the inequality 
OO) 


» 


Jyk 


ct (3.15) 














Ee, 


holds, where the constant C, depends on M. 
It is evident that by virtue of (3.2) for any point x, € 2, there exists a sphere 
of radius C,h with center at the point x, lying completely inside the surface 2, 


oe an 
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(C, depends on M). By the maximum principle the values of w (x) at the 
boundary of this sphere do not exceed one in modulus. The assertion of the 
Lemma follows from the above and from known estimates for derivatives of 
harmonic functions at the center of a sphere of given radius in terms of its 
values at the boundary of the sphere. 


Lemma 3. Let w (x) be a harmonic function regular inside the surface 2, 
of Lemma 2, satisfying the following conditions. 


1) For some vector x, € 2, and for an arbitrary point x lying on the surface 
2, and inside the sphere 


Ix—Xol SA, 
the function w (x) is 
o(x)=1. 
2) For an arbitrary xe 2, 
|w(x)|<1. 


Then the following inequality holds 
h \@ 
1—@(Xo,) <C, (4) (3.16) 


where C; are some constants depending on the constant M; x9, — is the 
vector with the polar coordinates 


Co—h, Go, Oo 


where 0, Yo, 99 are the polar coordinates of the vector xo. 

From the conditions imposed on the function f(@, 9) and also from the 
maximum principle it follows that the assertion of the Lemma is equivalent to 
the following assertion. 

Let W’ (x) be a harmonic function regular in the sphere with a cut-out cone 


pe: F-02e, 
(0, @, 9 are polar coordinates of x — 7 and & is a constant depending on M). 
Suppose W’ (x) takes on the following boundary values 
W'(x)=-1 
on the surface of the sphere e = A; 
W'(x)=1 


on the surface of the cone z/2 — @ = «. Then at the point x,, where the polar 
coordinates of the vector x, — 7 are equal to 
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Te 
Q,=h, o,=9, =z 


the function W’ (x) satisfies the inequality 


h\? 
1-wnse,(4) 5 


where C; are constants depending on «. 
The latter assertion follows from known estimates for harmonic functions. 


Lemma 4. Let w (f) be a nonincreasing function of the parameter ¢ on the 
interval [0, C] 


o(0)=1, w(C)=0. (3.17) 


Then for an arbitrary t,0< t<1,n > 1 thereexists anumber4,1< k<n 
such that 


o()—a(t)s—. (3.18) 


We consider the sum 


BE [o(t*)—o(t]J=0(")—o(0). 


From the conditions of the lemma all the terms of the sum are non-negative 
and the whole sum does not exceed unity, from which it follows that some 
member of the sum does not exceed I/n, i.e., for some & (3.18) holds. 

We now turn to the proof of the theorem. We introduce the following 
notations. 

Let 


f(g, = mar [F;(9.9)], 
fi(@, 6) = zat LS; (Q, 6)] : 


We denote by 2? the surface defined by the equation 


o=f'(9,8) 
and by 2! the surface defined by the equation 
o=f'(9,8); 


and by x ny the surface defined by the equation 
o=f' (9,0) +h. 
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We denote by zh the part of the surface a) whose points satisfy the inc- 
quality 


f1 (9,9) = fr (9, 9); 


and by ain we denote the part of xn) whose points satisfy 


fi (9, DS fr (@, 0). 


The corresponding parts of the surfaces £?, 5% will be denoted by 
pit gle gis yin, 
We denote by bah , the part of pe whose distance from a is not less than / 
r(x, LG) =A KEL is 


we denote by PD (h) A the part of & oid which belongs neither to Dares y nor to 
ma , and finally, by pa ay a2 We denote the part of Dar y Which does not belong 
to the surfaces 
fe f 
Sapa xii ae Zcnya, (42 >A,)- 


The parts of X*, 2? to which correspond rays from the point 7 intersecting 
the indicated surfaces, will be denoted by 


Pha = ie Ded A2 
De ke ee ees 
respectively. 
We denote by 2 the domain enclosed between the surfaces 5, S* while 
by QF, QF we denote the parts of 2 enclosed respectively between the surfaces 
af ae and pe pie and the radial straight lines bounding them. 


Finally we denote by u (2) the volume of the domain 2, by o (2) the area 
of the surface 2, by « (x) the function defined by the relations 


a(x)=1 xEQ5; 
a(x)=—1 xED, 5 
a(x)=0 outside Q; 


and by u(x) we denote the difference 


u(x)=u, (x)—1t2 (x). 
Clearly 
ao) 4 
r(x, as 


Let v (x) be a harmonic function regular inside the surface 2? and contin- 
uously differentiable on 2*. Then from potential theory 





u(x)= f 
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Ja(x): ooparaze j(2 urv— 2S o-u)de (3.19) 


where 0/0n is the derivative along the normal to 2. 
We set 


»(x)=div[W (x) («—n)] 
where W (x) is a harmonic function regular inside pai and assuming on on) 
the following boundary values 


W(x)=1 xe GA; 
W(x)=-1 xEZGy4; 
W(x)=0 xELiy4 


We estimate first the right-hand side of (3.19). By virtue of Lemma 1 and 
the conditions of the theorem, the following inequality holds on the surface 2% 


|u (x)| + |grad u(x)|<e,=C, |In Sica (3.20) 
and by virtue of Lemma 2 
|v(x)| + |grad v(x)|< a (3.21) 
From (3.20), (3.21) and (3.2) we obtain 


0 0 
J | ou »—So-udo 


where Cy is a constant dependent on M. 
We transform the left-hand side of (3.19) as follows 


f a(x) v(x) dx= f a(x) div[W (x)(x—n)] dx 


= J a-[W-(x—n)],do— § o- [W-(x—n)],do 





Cc C 
Ss eo (2) S75 (3.22) 





a ae 
+ J a[W-(x-n)],do— jf a-[W-(x—n)],do 
as ra 
+ J a [W-(x—n)],do— J a: [W-(x—n)],do 
per ae 
+ J a-[W-(x—m)]do— J a [W-(x—n)],do. (3.23) 
xi DH 


where [W. (x — 7)], is the component of the vector W (x). (x — y) normal 
to Landi > 4h. 
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Let us divide the terms of the right-hand side of (3.23) into three groups 
and estimate each group separately. 

Let us consider the first group. By virtue of Lemma 3 the function W (x) 
on the surface a, Pe differs from +1, —1 respectively, by no more than 


h \° 
Cs (4) ; 
A 
where C; are constants dependent on M. Therefore 


J a [W-(x=n)],do— J a-[W-(x—-n)],dot+ J «-[W-(x—-n)],do 


Dy a ig 
— J «[W-(x—n)],do 
pea 
7 h \&s 
> w(Q7)+ u(Qz)—Cy (+) : (3.24) 


It follows from (3.2) and the definition of Da , that if the straight line 
defined by the angles ¢, 6 intersects pa , then the following inequality holds 


Ifs (9,0) —f2 (9, MIS CA. (3.25) 
It clearly follows from (3.25) that 
1(Q)— p(Q7)—B(Q7)SCo°A (3.26) 
where Cy depends on M. Substituting (3.26) in (3.24) we obtain 
J a-[W-(x—-n)],do—- J «- [W-(x—n)],do 





pa bg 
+ fa-[W-(x-mhdo— § «-[W-(x—n)],do 
Oa bag 
h \s 

>u(Q)—C,A-C, (4) ; (3.27) 

We consider the second group. Clearly 
J a-[W-(x—n)],do— J a-[W-(x—n)],do <a (2; ,) +0 (2h) : 
Dy a Zh a 
(3.28 


We consider the function 


wo 2 IFA E=9(E) 


a(3')+0(2') 
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The function @ (A) satisfies the conditions of Lemma 4, and thus for any 
¢ < 1 and n there will be a number & such that 


w(t) —ot<—, k<n. (3.29) 


(<4) (3.30) 


Substituting (3.30) and (3.29) in (3.28) we obtain 


We set 


J oLW-(x—m]do— f a-[W-(x—n)] do 


eC gs (3.31) 
Lha Zh 





Finally we consider the third group. It follows from Lemma 3 that the 
function W (x) on the surface 2/, 3? satisfies the inequality 


h \Ow c 
IW(x)lsCry Tt =C,,t", (3.32) 
from which it follows that 
J «-[W-(x—n)],do— J a-[W-(x—-n)],do<C,3t%. (3.33) 
zi xi 
Substituting (3.27), (3.31), (3.32) in (3.23) and substituting (3.22) and (3.23) 


in (3.19) we obtain 
C48, 


2k+2° 
t 





C 
H(Q)S Cyt" + Cot! +— + Carts (3.34) 
In (3.34) we set 

t=— (3.35) 


where 7 satisfies the relation 


1 \2"*2 1 2n-1 
(ty sas{ ty. (3.36) 


Substituting (3.35) and (3.36) in (3.34) we obtain 





C 
1(Q) ST 


from which, by virtue of (3.2), the inequality (3.3) follows. 
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§ 2. A Class of Nonlinear Integral Equations 


Let us consider the equation 
Ag=f (3.37) 
where @, f are continuous functions of x on the interval [0, 1] and the operator 


A can be represented in the form 


1 g() 
Aga J P(x,¢,m)dndé (3.38) 


with P (@, &, 7) continuous and non-negative. The operator (3.38) is a special 
case of the URYSOHN operator 


Ag= PK (x,¢,9(6))d¢ 


when the kernel K (x, &, 7) is a continuously differentiable function, satisfying 
the condition 


K(x, ¢,0)=0 


Let us formulate a uniqueness theorem for equation (3.37) for a class of 
operators A having the representation (3.38). 


Theorem. Let the operator A satisfy the following conditions. 


1. For any 1,02, Q2 (x) => @, (x) satisfying a LipscuiTz condition with 
some constant C, the operator a Ay, exists where A, is the FRECHET deri- 
vative of the operator A at the point g. Moreover, for any w satisfying a Lir- 
SCHITZ condition with the constant C, the function A>) A, W satisfies a Lir- 
SCHITZ condition with the constant C, depending on C. 


2. If pz (x) > @, (x) the operator ee om can be represented as 


Aj, AgW = 5 Qo, [x, ae P2 (2)-¢ (€)] w(é) dé 


where the function Qo, (x, &, 7) is non-negative and continuous for 7 > 0. 
3. The integral 


iil 
J f Q, (x, 6,0) dodn 


is convergent for any and x. 


°o 


4. There exists a continuous function of two variables w (t, h) dependent 
only on the constant C and satisfying the conditions, 
a) w(t, 4) 9 0, hk > 0, t+ O; 
b) for any @ satisfying a LrpscHiTz condition with the constant C the ine- 
quality 
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hO(xEh) 
“TG hy SOCEM 


holds where 
Axth 


Wel h)= J J Opa Eh) dgdn. 


Then the solution of equation (3.37) is unique within the class of functions 
satisfying a LipSCHITz condition. 

It should be noted that this type of conditions is satisfied by the integral 
operators arising in the inverse problems of potential theory. 

The classical example of an inverse problem is the one arising in the theory 
of NEWTONIAN potential formulated by Novikov. Solving the planar inverse 
problem of Novikov is equivalent to solving equation (3.37) with 


P(x, €,h)=In[r? +h? —2rycos2a(x—£)-n], O<o<r. 


If p, (x), G2 (x), and W (x) are continuous functions, @, (x) > @, (x), then 
in this case the function A; Ay .W is a density of a single layer potential with 
weight @, (x) distributed on the curve defined by the polar equation 


= 2 (x) (2) 
(2 xx is the polar angle) and satisfying the following condition, the logarithmic 
potential of a single layer with density g, (x) & (x) on the curve 

e=9,(x) Wy) 


coincides with the single layer potential with density AS A,,w outside the 
domain bounded by I’. 

It is easy to show that in this case the operator Ay) A 4, satisfies the condi- 
tions 2, 3, 4 of the theorem. 

Let us introduce some notation. Let g; be continuous functions and let 


u(X, Py, @2)=max —,(x), 
j=1,2 
v(x, Oy, @2)=min —; (x), 
W (x, 1,92, 3)= min U(X, 9; O)= max D(X, @;, P;)- 
: 3 


i, f=1,2,3 i, j=1,2, 
Leta; (j = 1, ...) be sets of points of the interval [0, 1]. We denote by 0 (x1, %) 
the distance between the sets «,, x». 
We denote by « (@,,@,) the set of values of x for which 


91 (X)> @2 (x) 
and by & (@;,@,) the set of values of x for which 


P1 (x)=@2(x). 
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Let A; = 1, ...) be some real numbers. We denote by « (~,@2, 4) the set 
defined by the relations 


, xX€4(1,2,A) 
1 


xEa(P, (2); o[x,a(91, @2)| =A 
and by &(@1,Q>, 2) the set 


a(P,, Pa> 4) =a(,, P2)OK (GY, 92)0«(G,, Pa, A) 
while by &((;,@>, 4,, 4) the set 


&(P1, G2, 41,42) =4 (Py, 2,4, OX (Py, Pz, 42). 


The proof of the theorem is based on the following lemma. 


Lemma. Let A be an operator satisfying the conditions of the theorem and 
let p, and @, be functions satisfying a LrpscuiTz condition with the constant C, 
and such that the sets « (~,,@,) are not empty. Then for any sufficiently small 
positive A and / there exist functions @,, @,, W possessing the following prop- 
erties, 


L. u(x, G1, G2) 2G; (X)20(%, 91,92), J=l2; 
2. Ag,—AQ,=AG,—AG, + AW, U=U(X, G1, G2); 
3. W(x)20, xEea(P2, 9); 
w(x)s0, xEa(P,,92); 
4. e[4(1,92), «G2, G,)]=4; 
max | (x)— (x)|=h; 
6. the functions W, @; (j = 1, 2) satisfy a Lipscuitz condition with the 
constant C, depending on C. 


We prove the above lemma. Let ~,,@, be functions satisfying a LrPscHITZ 
condition with the constant C, let 


Ag, —Ag,=f 


and let 2, A be some sufficiently small constants. 
We denote by 7, the function defined by the relations 


1 = 
a) Vmi== (2-91): xE&(P2, 91,4); 


1 eo 
b) Ymar= (1 — 2); xE4(P1, 0254); 


Yma=O, XEX(Py, P2,2A)\Ga(G2, 91,24); 
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¢) Y ma iS a continuous function on the interval [0, 1) and it changes linearly 
on the set &(@,,@p, 24, A). 
We consider the sequences of a function Q mza>P2mka Wma (K = 1, «Dy 
P jmoa = Pj } (j=1,2) 
P jmk+12=Pjmkat 9 jmka 


Wmoa=9, Win, kt 1,2= Wma t mea» 
where 


A) Ota =2 Omar Pamka— Pimka)s XEa(P2, Px); 
5.2mka=P mas Pimka— Pamka)s X€0(P1; P2); 

b) Oma = W [ARS Al a0 en 0; Prmca—Pimeals X€4(P1, Pr); 
S2mka=W [Abi Avie meas 0; Pinka—Pamkal, X€a(P2,91); 
Umja (X)= U(X, Prmkas Pama) > 


Umka (x) =v G; Pimka> Pomna) > 


Smnzx (X) = 5 mea (%)> XEa(P2, 1), 
Bnza(X) = Samia (X), X€4(Q 1, 92); 
C) 0 jmza=9, XEX(P1, P2); 
d) Senta =A yarns mba — Ot mas xEa(Q1, M2), 
Oink = Abus vmna® mka — Samar X€E4(P2, 91), 
Srna = Ausra A vmsca © mkh? XEX(P1, 2). 


One can easily see that 


C,k {1 
40am A nts Aantntr SIs (1) (3.39) 


where C, is a constant depending on C and yz (£) is the modulus of continuity 
of the function P (x, é, 7). 
Indeed by virtue of the definition of the FRECHET derivative 


AP jmk+1a— AP jmka— A 6 jnza imkal =0 (119 jmaall) . (3.40) 


From the representation (3.38) of the operator A it follows that the right- 
hand side of (3.40) satisfies the inequality 


C 1 
0(I5 pals x( 4). 


m 


The inequality (3.39) now follows by virtue of the definitions of the functions 
Pj mk Oj mkas W mka: 
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Next, from the conditions imposed on the operator Avi Ag, it follows that 


for k < m the functions Qj mya,Wmxa Satisfy a Lipscuirz condition with a 
constant C3 dependent only on C and that fork > m 


P jm, k+1,2= Pjmka=Pjma>s 


Win, c+ 1,2=Wmka= Wma - Gay) 
The functions @;,,, in (3.41) evidently satisfy the conditions 

21% (Pima Pama) “(Pama Pima ZA, 

0 (91, 2) SP jma SU (G1, G2) oy) 
and for any sufficiently small A, A 

max |@2ma—Pimal >2h. (3.43) 


We consider now the sequences of functions ; ma,P2mar WU ma (mM = 1,...). 
By virtue of the fact that all the elements of sequences are bounded and satisfy 
a Lipscuitz condition with the constant C, independent of m, one may choose 
from the sequences uniformly convergent subsequences 

Pimgas Pamga Winga (q=1, Ee 

By virtue of (3.42) and (3.43) the functions 


G;= iM @jngas (j=1,2) 


qr a 


"7 = lim Waa 


qr on 


satisfy the following relations 


o[«(G,, G2), a(G2,G,)|24, 

AG,—AG@, +AW=f, 

v(P1, P2) Gj SU (G1, G2), (3.44) 
max|@2— ilo" 


and a Lipscurrz condition with constant C4. 
We denote by y,,, the function defined by the relations 


ae a0 
¥m=— (G2—P1)s XE4(G2, G1), 


Ym=—(G1-— 2), x€a(G1, G2) 


Ym=0 x€&(G1, G2) 
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and consider the sequence of functions Qj mx; P2mko W mk (K = 1,...) 


Q jmo= Pj (j=1,2), 
P jm, k+ 1 =P jmk FO jmk » 
Wino =W 3 


Wm, += Wink + Ont 


where 0; mx 5 mz are functions defined by relations similar to relations a), b), c) 
defining the functions 6; mas 5 mka- 
In complete analogy with (3.39) we obtain 


Czk 1 
|A@am~—A tmk+ Aum mk — J | =o, (1), (3.45) 


max |P2mp— Pimpl >h, 
MaX |Prmp+1—Pimp+il Sh. 


By virtue of the definition of the sequences 9;,,, we can obviously find a 
number p, 0 < p < m, satisfying (3.46). 

We denote by G;,,, W the functions Djm = Pjmp» Vim = W mp and consider 
the sequences 


(3.46) 


D jms Vm (m=1,...). 


In analogy with the sequences 9,,,;, Wing (7 = 1,...) we may select uni- 
formly convergent sub-sequences Pj m,? WV mg (q = 1,...) from the sequences 
Dim Wm: By virtue of the above, the functions 

Q;= lim P jmg (j=1,2), 
qr a 
w= lim ¥,,, 
q7 no 
satisfy the relations 1)—6) and the lemma is proved. 

We now prove the theorem. We suppose the solution of (3.37) to be non- 
unique within the class of Lipscuirz continuous functions, i.e., there exist 
distinct functions @,,@, satisfying a Lipscuirz condition with the constant C 
and 


Ag,—Ag,=0. (3.47) 


It follows from the Lemma that for any sufficiently small / and 4 there exist 
functions @,, Mp, W satisfying the conditions 1)—6) of the Lemma and there- 
fore, according to (3.47), 


AG,— AG, +A, W=0. 


Applying the operator Az} to (3.48) and making use of the conditions 2) 
and 3) of the theorem, we obtain 
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1 
J J OxLx,¢.n] dn-sign[9,(6)—G, ()}dE+P(x)=0. (3.49) 
We denote by X a value of the independent variable x such that 


|P2(X)— G; (X)l=max |G. (x)— , (x)]. 


Because of the properties of @; 


G2 (X)—@, (X)=h. (3.50) 
Let, for the sake of definiteness, 
P2(X)>@, (X), 
1.€., 

P2(X)— G(X). (3.51) 


We put x = X in (3.49) and divide the integral in the left-hand side of (3.49) 
into three parts 
X+hy G(G)—v (G) : 
Q;(X, €,) dy sign[ Oz (C)—@, (Q] dé 
<= VY 0 


X—hy u(S)—v (6) 


+f J Q,(X, ¢, h) dy sign |@2— 1] dé (3.52) 


0 
1 a@(3)—e (8) : 
“F j j 0; (x, G3 h) dy-sign [9, —9,] dé+W(x)=0, 
xXx+ 10) 


Ay 
where 


h,= c 
ae (ore 

We denote the integrals on the left side of (3.52) by J,, J, and J3 respectively 

and estimate them separately. Let us consider the first integral. From (3.51) and 

condition 6) of the lemma it follows that 


7 - h oe, 
lP2)- PL O25 Ix— X27 om (3.53) 





from which 
I, >W,(s,h), 


where W, (x, h) is the function defined in condition 4) of the theorem. 
We consider now the second and third integrals. From condition 4) of the 
lemma we have, 


X¥-Aa(S)—v (8) 


I> ff Q,(%,é.n)dn-sign[G,-G,] dé, (3.54) 
0 0 
1 a(€)—v(g) ; _ 7 : 
=f § Q,(%,é.n)dn-sign[¢,—9,]4é, 
X+A 0 


5 Springer Tracts, Vol. 1!, Lavrentiev 
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while from (3.54), condition 4) of the theorem and condition 5) of the lemma 
we obtain 


I, +13> —@(A,h)- W,(x, hh). (3.55) 


From condition 3) of the lemma and from (3.51) it follows that the last term 
in the left-hand side of (3.52) is non-negative, i.e., 


W(x)>0. 


From (3.54), (3.55) and (3.56) it follows that for sufficiently small # there holds 
the inequality 
T,+1,+1,+W(x)>0 


contradicting the equality (3.52). Thus the theorem is proved. 


§ 3. Inverse Problems for Some Non-Newtonian Potentials 


Inverse problems for NEWTONIAN potentials have unique solutions only 
under very essential restrictions on the nature of the density distribution. An 
example constructed by P. S. Novikov indicates that the solution of the inverse 
problem of NEWTONIAN potential may be non-unique even in the class of 
densities which are piecewise constant with a given constant and are contained 
in simply connected domains. An analogous situation seems to take place 
also for potentials corresponding to the fundamental solutions of other elliptic 
equations. 

As will be shown in this section, for inverse problems with some Non-NEw- 
TONIAN potentials not connected with differential equations, uniqueness holds 
for much milder restrictions on the nature of the density; its continuity and 
finiteness will be quite sufficient. 

Let x, € be n-dimensional vectors with components (x, ..., Xp), (&4) «5 Fn)» 
let D, D, be non-intersecting simply connected bounded domains of n-dimen- 
sional space and let 





1 
Oe ren OO (3.56) 
where 
r(x, €)=|x—€]. 


Given the function u(x) in D,, it is required to reconstruct the potential 
density o in D. 


Theorem. For p > n + 2 the solution of the problem is unique in the class 
of continuous ¢ (x), i.e., the continuous function o (x) is uniquely determined 
by the values of the function u (x) in the domain D,. 
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For our proof we shall examine the function 


_ 
v(x)= ee mee epee (3.57) 


R(x, O= Jr, H+ ly" 


and y is a vector with components ()',, ..., ¥m+;), m= p— un — 3. 

The function v (x, y) will clearly be a harmonic function of the variables 
x, y with cylindrical symmetry. From the results of § 6, Chapter II it follows 
that the function v (x, y) is uniquely determined on a set Q of x, y space: 


where 


x, yEQ, if xXED or y +0; 
by its values on the set 
xED,, y=0. 


The function v(x, y) is a NEWTONIAN potential of a simple layer with 
density @ (x) distributed on the n-dimensional manifold x é D, y = 0 in the 
(1 + m+ 1)-dimensional x, y space. It is clear that 9 (x) is uniquely deter- 
mined from v (x, y). Hence the theorem is proved. 

We note that the method used for the proof of the uniqueness theorem 
permits one to obtain estimates characterizing thc stability of the solution. 


§ 4. An Inverse Problem for the Wave Equation 


In this section we shall prove a theorem giving the uniqueness for the 
solution of an inverse problem for the wave equation. 

We consider the equation 

2 
= du (3.58) 
where « is a function of the three variables x, y, ¢ and # is a function of the 
variables x, y. 

We consider the following problem: the domain Dg is given in the x, y 
plane. The function n(x, y) > 0 is continuous and identically equal to one 
outside of Dj. Moreover, in some domain D, with D, M Dy empty, there is 
given a family G of solutions of (3.58) for all ¢ > 0. It is required to determine 
the function 7 (x, y) inside Dp. 

We note that a similar problem in the case of one variable was considered 
in the paper [59]; quite close to the problem being considered there, is the well- 
known STURM-LIOUVILLE inverse problem. 

We formulate a uniqueness theorem for the above problem for a certain 
family G. 
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Theorem. Let D, and D, be bounded and simply connected domains, and 
D,,a bounded, simply connected domain, which does not intersect with Dy or 
D,. Let G be a set of solutions of (3.58) satisfying the following initial condi- 
tions 


u(x, y,0)=0; 
0 
ar” (x, y,0)=6(x — x9, Y— Yo) (3.59) 


where Q (xy, Yo) is any point of Dz. 

Then the solution of the inverse problem is unique, i.e., the function 
n (x, y) is uniquely determined within Dp. 

We give a brief sketch of the proof. We denote the solution of the CAUCHY 
problem (3.59) for equation (3.58) by u (x, ¥, Xo, Yo, 2) and consider the function 


DIX Vi Xos Vo R)S= f U(X, Y,Xo, Vo, t)cosAtdt. 
One can easily see that ; 
Av= —8(x—X9, y— Yo) —A?n7v. (3.60) 
The function v is the fundamental solution of the HELMHOLTZ equation 
Av=—2?n?v 


with the singularity at the point Q (Xp, yo). It is well-known that a fundamental 
solution for an elliptic equation with analytic coefficients is an analytic function 
of both the independent variables x and y and the coordinates of the singularity 
Xo, Yo everywhere away from the singularity. 

According to the conditions of the theorem the function v is given in the 
domain D,, of the 4-dimensional space R (x, y, Xo, ¥o), the direct product of 
the domains D,, D, of the spaces P (x, y) and Q (Xo, yo). The function 7 (x, y) 
is identically equal to one outside of Dj. Hence, by virtue of the uniqueness of 
the analytical continuation, the function v may be considered to be given 
everywhere outside of the domain Dog = Dy X Do. 

Let D3 be a bounded domain containing the domains Dp, D, and D, and let 


a(x t(yon)’, ro=(%o— 6) + (Yo —0)’. 
From (3.60) it follows that for Pe Dg 


1 
v(x, Y,Xos Vos A)=5— In [(x—x)*+(v—yo)7] 


—5—An? (x, y) J 0m Xor Yo°A)In r de dy 
Tt Do 


+5(X, Y, Xs Yo. €)3 (3.61) 
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a oe 3 Ov 
are (2m } —5, In r)as, 


where /’3 is the boundary of Dg. 
We denote by 2, (x, y, x, ¥p) the function 


6) bee Ys Xo. Vos ee 0° (X, Pe tas Vos | 





0, = 
: OA Ox: OXo Oy . Oyo 


Tt can be easily shown that the function v,, for RE Dg3 = Ds Ds, equals 


v= [mG DOD ae ray ss, (3.62 


where 





eas ee 076 - 075 
Lye) POET OL = GA. Oecaxa, -OysOye | 
By virtue of (3.62) the functions v,, 5, for RE Dgg and Re Dog are analytic 
functions of the variables x, y, Xo, Yo- 
By setting 
X=No, Y=Jo; P(x, y)eD, 
in (3.62) we obtain 


s 1 7 
v2 =(0, (x,y, x, y)—8, (x, », x, y)= J n (En) —dedy . (3.63) 
Do 


The assertion of the theorem follows immediately from the results of § 3 
and (3.63). 

We note that the theorem generalizes to the case of the wave equation in a 
space of an arbitrary number of dimensions as well as to the heat equation and 
some equations of hyperbolic and elliptic type of higher orders. 
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In this section a statement more general than that of the previous section 
is considered. A class of inverse problems is reduced to linear integral equations 
of the first kind which are then investigated by methods differing from those 
of the previous section. 

We consider the equation 


0 0 
Pe (2) u(x, y)=P3 (2) Ps (2) u(x, y), (3.64) 


where x, y are vectors with components (x;,...; Xn), 4, ++. Yq) and P; are 
polynomials with coefficients continuously depending on x. 
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In regard to P,, we assume that the equation 


(2 Joa=0 


has a fundamental solution G (x, x°) defined in the entire x, x) space*. 

We pose the inverse problem for equation (3.64) which will be considered 
in this section. 

Let Dy, D, be some bounded domains of x space, where the intersection 
D, © Dg is empty. The coefficients of the polynomials P,, P, are given for the 
entire x space and the coefficients of the polynomial Ps are given everywhere 
outside of Do. In addition, let a family of solutions u (x, y, €) of (3.64) depending 
on the parameter & (&,,..., &,), be given in the domain D,. 

The functions u (x, y, 5) satisfy the conditions 


r 


——— u(x, 0, =f, --+ x, (% 6), 
Oye... dyke ( y= Si : 

0"'u (x, 0, ¢) 
Oxi! Ox! 


Ou (x, y, €) 


kn 


Oyt! ... Ovi 


where the indices k,, ..., k,, run over the values which are less than the maxi- 
mum values of the corresponding indices in the polynomial 


0 
= (2) 


and i, + 4,,-..,/, + J, run over a similar set with respect to the polynomial 


0 
Ps (2). 


ee is a sphere of radius 0 with the center at the point af (x, &) are given 
finite functions, and C,, C, are some constants. 

It is required to determine the coefficients of the polynomial P3 in the 
domain Dg. 

In the case when 7 = i = 1 the above inverse problem is equivalent to 
the known SturM-LIOUVILLE inverse problem for ordinary differential equa- 


tions. 


3G (x, x°) d 


Oxi! on ots oa oe) 


lim j 


27a Sexo 














<C,+lyl? 





* The existence of a fundamental solution has been proved, say, for elliptic equations 
of the second order and also for ordinary regular equations with constant co- 
efficients. 
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We denote by v (x, A, &) the function 


oO 


v(x,4,6)= J u(x, y, dyexp{—(A, y)} dy 


0 
(A is a vector with components A,, ..., 4,,)- 


By virtue of (3.64) and (3.65) the function v (x, A, £) satisfies the following 
differential equation 


3 
Pi(2)omPs(—a)-Pa(2 oa, (3.66) 


where J 


oe) 
W(x, A, €)=Q (2 2) u (x, 0, ¢) 


is a function which appears as a result of integrating by parts the expression 


for 
Ce) 
P,{ — }- 
(2) : 
and Q is a polynomial. 


We multiply both sides of (3.66) by the function G (x, x°) and integrate 
over the whole x space. From inequalities (3.65) and (3.66) it follows that 
the integrals will converge, and by integrating the corresponding expressions 
by parts, the boundary terms at infinity will be equal to zero. 

Thus, we find that the function v (x°, 2, ) satisfies the integral relation 


, 6) 
v(x°,4,2)= | G(x, x°) P3(4,) P2 (2) vdx+ J G(x,x°) pdx. (3.67) 
Jj 
From (3.67) it follows that the function v (x, A, €) is an analytic function 
of A in some neighborhood of the origin. We consider 


r 


Ce) 
Ux, th ky (x, aan age 0, é), 


OD csOAe 
TR Co eee 
BPD SEME NOS Ortega eames es ot 
Oo 
P,... P,(0). 


re ake. Odkm 


“mm 
By virtue of (3.67) the functions vy, --. go> Why e+ & 
following system of relations 


ie OO) = [CC “| See ree (2). Uf ax coh dx 


+ f G(x, KV, eeu XyC) aes (3.68) 
vy. 9X O= J G(x, x) Ho. oa, Adz. 


Pia’, 


1 


. k, Satisfy the 


im 


68 Inverse Problems for Differential Equations 


Relations (3.68) may be regarded as a recursive system of linear integral 
equations of the first kind in D, for the coefficients P;, ... ; of the polynomial 
P3. In fact it follows from (3.68) that the function ug... eles &) may be re- 
garded as given in the whole space. 

We consider now the function 


6) 
Vo. age: = Ge, x°) Po .1...0P2| a— } 0... 0% 6) dx 
—— Do —“~—— Ox; 
q q 


re) 
+ J G(x, x°) Po 230Pa( ee) v6 no 0d 


Do ¥ 


+ J G(x, x°) Wo ..1...0(% €)dx. (3.69) 

eo 
If x®e D, then in the equality (3.69) all the terms may be regarded as given 
with the exception of the first one in the right-hand side, and therefore (3.69) 
may be regarded as a linear integral equation of the first kind for the function 


Praia): xEDo 
“~— 


q 
with the kernel 


fe) 
K(x, x°, )=G(x, x°) P, (2p. 0 (% €), x°eD,. 
Jj 

If the above equation has a unique solution for any qg then the equations 

defining the function 
Po te! 1 Oe 1-9 (X) 
—_" 
q a 

is found in a similar way, and so on. 

We will now carry out a detailed consideration of the system (3.68) in the 
case in which the original differential equation is of the form 


B J 
Aiu(e y=). a, Oren (x,y) (3.70) 


( is a scalar, «, 8 are integers.) 
The system (3.68) in the case under consideration will take the form 


5, (x°, €)= J G(x, x°) v9 (x, €) ay (x) dx, 
09 (x°,€)= i G(x, x°) Wo (x, é) dx, 
5, (x°, Y= ane J G(x, x°) oj (x, &) a;(x) dx +0, (x°, 8 
a G(x, x°) ay (x) 09 (x, €)dx— J G(x, x°)- W(x, O dx, 


(3.71) 
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Yan LX —x°|77 "In |x—x?], for n even, 2a—n>0; 
G(x, x°)= 


Yen bax 


in the remaining case. 

From (3.71) it follows that for uniqueness of the solution of the inverse 
problem under consideration it is sufficient that there be a unique solution of 
the integral equation of the first kind with respect to the function a (x) 


J G(x, x°)-v9 (x, &) a(x) dx =a(x°, &). (3.72) 


We multiply both sides of (3.72) by an arbitrary function p (x°) and integrate 
over the domain D,. 


[O° DeG)dx?= [ W(x) v0( Aa(x)dx, 


W(x)= J G(x, x) 0 (x) dx?. ved 
Dy 

The function W (x) in (3.73) is a potential with density 0 (xg) distributed 
in the domain D, whereas the function v, (x, £) is a potential with the density 
Wo (x, §) distributed in the domain D,. Hence, for uniqueness of a solution 
of (3.72) it is sufficient that the linear hull of products of these potentials should 
be dense in the set of functions defined in Do. 

Let, for the moment, the set of functions Wo (x, &) be such that the set of 
potentials vy (x, €) is dense in the set of all harmonic functions regular in some 
extension of the domain Dy — Do, *- 

From potential theory it is known that the potentials W(x) lie densely in 
the set of all harmonic functions regular in Dp,. It can be easily seen that in 
this case the linear hull of defined products of potentials is dense in the set of 
continuous functions defined in Dy. In fact, we consider the harmonic poly- 
nomials of order v which are normalized and homogeneous with respect to 
the coordinates of the points x — x! where x! is an arbitrary point of Dp. 
Tt is well-known that the sum of squares of all these polynomials is equal 
to |x — x1|® whereas the function {1 — |x — x*|®/R?} approximates 6 (x — x1) 
in the sphere |x — x!| < R. Thus, we have proved the following uniqueness 
theorem for the inverse problem for equation (3.70). 


Theorem. If a family of solutions u(x, y, £) of (3.70) is such that the set 
of potentials v, (x, &) is dense in the set of all harmonic functions regular in 
the domain Dg, (some extension of the domain D,), then the solution of the 
inverse problem for (3.70) is unique in the class of continuous functions a, (x) 


(C24, 2468): 
* This condition is fulfilled if, for instance 
Yo (x, €) = u (x, 0, €) = 5d (x— é) 


where the parameter & runs over the domain D,. 
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